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NOMENCLATURE 


exponent  for  along-wind  variation  of  concentration  fluctuation 
intensity,  see  (3-4) 

crosswind  spread  exponent,  see  (A2) 

gas  molecular  weight 

3 3 

time  or  ensemble  averaged  mean  volume  fraction,  m /m 

3 

time  or  ensemble  averaged  mean  mass  concentration,  kg/m 

2 6 

concentration  fluctuation  variance,  kg  /m 

effective  isokinetic  source  diameter  of  a vertically  directed  jet, 
see  3-21,  m 

"source"  diameter  of  a jet  after  expansion  from  a high  pressure 
release,  see  (3-22),  m 

source  diameter  for  isokinetic  injection  in  the  x direction 

crosswind  spread  constant,  see  (A2)  and  (A8) 

mass  fraction  of  pollutant,  kg/kg 

4 3 

buoyancy  flux,  see  (D4),  m /s 

4 2 

momentum  flux,  see  (D3),  m /s 

(h^  + Ah)  source  height  above  ground,  m 

thickness  of  the  lower  portion  of  the  atmospheric  boundary  layer 
that  is  influenced  by  mechanical  turbulence  from  surface  roughness, 
m 

plume  rise,  m 

buoyancy  rise  of  a plume  at  a specified  downwind  distance,  m 
buoyancy-induced  component  of  final  plume  rise,  m 
final  rise  of  a purely  buoyant  plume,  m 

actual  final  rise  of  a plume  including  buoyancy  and  momentum 
components,  m 

approximate  final  rise  estimated  by  combining  the  separate  rises  of 
a purely  buoyant  plume  and  a pure  momentum  jet  using  a sum  of  cubes 
law,  see  (5-12),  m 


momentum  rise  of  a bent-over  jet  a specified  downwind  distance,  m 


Ahpif  final  rise  of  a pure  momentum  jet,  m 


i 


i 

00 


i 


P 


k 


1 


k 


momentum-induced  component  of  final  plume  rise,  m 

effective  height  of  the  concentration  variance  source,  see  (3-14), 


total  RMS  concentration  fluctuation  intensity  in  homogeneous  flow 
(i.e.  in  the  absence  of  wind  shear  or  ground  surface  effects) 


RMS  concentration  fluctuation  intensity  in  the  absence  of  the  image 
sources  and  sinks  (which  account  for  dissipation  of  c'  by  wind 
shear  and  other  ground  surface  effects),  see  (3-13) 


fluctuation  in  ten  si  ties  formed  by  the  ratio  of  conditional  RMS  and 


? 0 5 

mean  values  * ^nd  c*p  which  include  only  periods  when  non-zer 


values  of  concentration  are  observed 


1)  wavenumber,  proportional  to  inverse  wave  length,  2Ti/wavel ength , 
m"^ 

2)  ratio  of  constant  pressure  to  constant  volume  gas  specific  heat 


integral  length  scale  of  atmospheric  turbulence  at  wavenumber  k,  m 


integral  length  scale  of  self-generated  plume  turbulence,  m 

integral  Lagrangian  length  scale  of  the  autocorrel  ation  of  crosswin 
velocity  v at  a point,  m 

integral  Lagrangian  length  scale  of  the  autocorrelation  of  vertical 
velocity  w at  a point,  m 

m mass  flow  rate  of  pollutant  and  entrained  ambient  air  through  a 

plume  or  jet  cross  section,  kg/s 


total  mass  flow  rate  from  a source;  including  pollutant  mass  flow 
rate  Q and  flow  rate  of  other  components  of  gas  mixture,  kg/s 

total  mass  released  from  a transient  source,  kg 

power  exponent  in  windspeed  profile,  U 

2 

absolute  pressure  of  the  atmosphere,  N/m 

o 

absolute  pressure  of  gas  at  the  pipe  exit,  N/m 
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steady  mass  release  rate  "source  strength"  of  pollutant,  kg/s 

plume  or  jet  radius,  m 

radius  of  pipe  exit,  see  Figure  4-1,  m 

radius  of  the  "source"  jet  after  its  expansion  to  atmospheric 
pressure  outside  the  pipe  exit,  see  Figure  4-1,  m 

time  after  start  of  release,  s 

weighted  release  duration,  see  (2-6),  s 

travel  time  to  actual  final  rise,  s 

travel  time  to  final  rise  of  a purely  buoyant  plume,  s 

travel  time  to  final  rise  of  a pure  momentum  jet,  s 

release  time  interval  over  which  mass  is  emitted,  s 

windspeed  at  height  z above  ground,  m/s 

root  mean  square  x direction  turbulence  velocity,  m/s 

average  convection  windspeed  through  a plume  y-z  cross  section,  m/s 

wind  speed  at  the  local  plume  height  h,  m/s 

RMS  turbulence  velocity  scale;  usually  u',  v',  or  w' , m/s 

root  mean  square  y direction  turbulence  velocity,  m/s 

entrainment  velocity 

root  mean  square  z direction  turbulence  velocity,  m/s 
vertical  mean  velocity  of  the  plume,  m/s 

mean  velocity  of  the  vertically  directed  jet  at  the  "source",  after 
pressure  outside  the  pipe  exit  adjusts  from  to  , m/s 

downwind  distance  from  source,  m 

t^U^,  release  length  of  a transient  plume,  m 

crosswind  distance  from  plume  axis,  m 

height  above  ground  level,  m 

reference  height  for  turbulence  scale  evaluation  in  a boundary 
layer,  see  (3-9),  m 


z . reference  height  for  turbulence  scale  evaluation  in  inertial 

subrange  model,  see  (3-9),  m 

roughness  height  parameter  in  log-law  velocity  profile,  m 

Greek 

a 1)  strength  of  the  concentrati on  fluctuation  variance  sink  used  t 

2 

model  surface  layer  c‘  dissipation,  see  (3-15) 

2)  al so:  entrainment  constant  for  a jet  with  no  crosswind 

bending-over,  see  Figure  4-1 

p entrainment  constant  for  a bent-over  plume  or  momentum  jet,  see 

Figure  4-1. 

2 3 

e dissipation  rate  of  atmospheric  turbulence  m /s 

_ 2 
e dissipation  rate  of  self-generated  internal  plume  turbulence,  m / 

<))2  0.706  0^9  vertical  displacement  of  the  effective  concentration 

fluctuation  variance  source  above  the  actual  source  height,  for  a 
ground  level  release 

y intermi ttency  factor;  fraction  of  time  that  non-zero  concentratio 

are  observed  at  a fixed  receptor 

intermi  ttency  factor  on  plume  axis  at  y = o and  z = h 

3 

p density  of  ambient  air,  kg/m 

P^  density  of  gas  at  the  "source"  condition,  at  atmospheric  pressure 

after  any  expansion  from  pipe  exit  pressure  has  occurred,  kg/m 

ensemble  averaged  instantaneous  crosswind  spread,  m 

ensemble  averaged  instantaneous  vertical  spread,  m 

meandering  component  of  total  crosswind  spread,  see  (2-1),  m 

meandering  component  of  total  vertical  spread,  see  (2-1),  m 

apparent  source  size  "spread",  m 

plume  spread  due  to  sel f-genera ted  turbulence  caused  by  jet 
momentum  and  buoyancy,  see  (4-11),  m 

total,  time  or  ensemble  averaged  along-wind  spread,  m 

total,  time  or  ensemble  averaged  crosswind  spread,  m 
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normalized  plume  spread; 


used  to  characterize  downwind 


integral  Eulerian  length  scale  of  along-wind  turbulence  velocity  u 
in  the  x direction,  m 

virtual  origin  displacement  used  to  account  for  the  effect  of  source 
size  on  concentration  fluctuation  intensity,  see  (3-1)  and  (3-2) 

virtual  origin  for  source  size  effects  on  concentration  fluctuation 
intensity  in  inertial  surange  model,  see  (3-4)  and  (3-5) 

virtual  origin  for  conditional  fluctuation  intensity,  see  (3-19)  and 
(3-20) 

(Av  + Ay^)/2,  average  integral  scale  of  turbulence  contributing  to 
plume  spread,  m 

integral  Eulerian  length  scale  of  crosswind  turbulence  velocity  v in 
the  X direction,  m 

integral  Eulerian  length  scale  of  vertical  turbulence  velocity  w in 
the  X direction,  m 

crosswind  plume  spread  caused  by  the  combined  effects  of 
self-generated  and  atmospheric  turbulence,  see  (4-12),  m 

vertical  plume  spread  caused  by  the  combined  effects  of 
self-generated  and  atmospheric  turbulence,  see  (4-12),  m 

total,  time  or  ensemble  averaged,  vertical  spread,  m 

decay  time  constant  for  a transient  release,  see  (2-8),  s 
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EXECUTIVE  SUMMARY 


This  report  is  the  fifth  in  a series  which  began  in  1979  with  "The  Release 
and  Dispersion  of  Gas  from  Pipeline  Ruptures".  These  reports  have  concentratec 
on  developing  analytical  methods  for  predicting  the  rate  of  release,  plume  rise 
and  atmospheric  dispersion  of  toxic  or  flammable  gas  from  a high  pressure  pipe- 
line rupture  or  well  blowout.  The  dispersion  model  which  has  been  formulated 
predicts  not  only  the  mean  concentrati on  but  also  the  probability  of  exceeding 
specified  peak  values.  This  report  deals  with  four  distinctly  different  aspect 
of  this  dispersion  model.  The  new  analytical  methods  developed  and  the 
conclusions  drawn  from  them  are  summarized  below. 

Puff  Versus  Plume  Models  for  Transient  Releases 

An  accidental  release  is  usually  a short  term  transient  event  that  lies 
somewhere  between  an  instantaneous  puff  and  a continuous  plume.  Puff  and  plume 
models  are  examined  to  show  that  during  the  early  stages  of  dispersion  the  long 
sausage  shaped  release  from  a suddenly  started  and  stopped  source  behaves  like 
plume.  Later,  as  along-wind  dispersion  entrains  atmospheric  air  into  the  head 
and  tail  of  this  plume,  the  release  behaves  like  a large  puff.  During  the 
transition  period  from  plume  to  puff  behavior  both  theoretical  models  over- 
predict the  concentration  in  the  transient  plume  and  provide  a conservative 
estimate  for  concentration.  Because  both  puff  and  plume  models  are  conserva- 
tive, it  is  recommended  that  the  predicted  concentration  be  taken  as  the  small- 
est of  the  values  computed  by  puff  and  plume  approximations  to  the  transient 
release.  This  approach  should  simplify  the  risk  analysis  procedure. 

The  Effect  of  Source  Size  on  Concentration  Fluctuations  in  Plumes 

The  diameter  of  the  source  plays  an  important  role  in  determining  the  peak 
values  of  fluctuating  concentrations  because  plumes  from  small  diameter  sources 
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tend  to  meander  back  and  fourth  in  large  scale  atmospheric  turbulence  and  pro- 
duce higher  peak  concentrations  than  those  from  large  sources.  The  most  impor- 
tant variable  is  the  ratio  of  the  source  diameter  to  the  scale  of  atmospheric 
turbulence.  A new  theoretical  model  for  the  effect  of  source  size  is  developed 
by  considering  the  interaction  of  the  plume  with  atmospheric  turbulence. 

Because  this  new  approach  makes  direct  use  of  atmospheric  turbulence  scales,  it 
should  have  more  general  applications  than  the  previous  model  of  Wilson  (1982) 
which  is  based  on  experimental  data  from  wind  tunnel  simulations. 

For  high  speed  jets  from  pipeline  ruptures,  the  emerging  jet  entrains 
ambient  air  and  expands  in  size  before  it  has  been  slowed  down  enough  to  be 
influenced  by  atmospheric  turbulence.  A simple  theoretical  model  is  developed 
to  deal  with  this  effect,  and  it  is  shown  that  the  effective  source  size  from  a 
high  pressure  rupture  may  be  a hundred  times  or  more  larger  than  the  pipe  exit 
diameter.  This  larger  effective  source  size  reduces  the  probability  of  being 
exposed  to  high  peak  concentrations  in  the  plume  far  downwind. 

Dilution  by  Sel f-Generated  Turbulence  Near  the  Source 

The  high  velocity  buoyant  jet  from  a pipeline  rupture  will  generate  turbu- 
lence due  to  its  momentum  and  buoyancy.  Air  entrainment  from  this  self- 
generated turbulence  causes  the  plume  to  be  diluted  even  in  the  absence  of 
atmospheric  turbulence.  Estimates  are  made  for  the  rate  at  which  this  dilution 
occurs  in  the  high  speed  jet  near  the  source,  and  in  the  plume  after  it  has  been 
bent  over  by  the  wind.  The  theory  predicts  that  as  the  plume  enters  its  bent- 
over  phase,  ambient  air  entrainment  by  self-induced  jet  turbulence  causes  dilu- 
tion of  a hundred-fold  or  more  in  source  concentration.  However,  atmospheric 
turbulence  eventually  becomes  the  dominant  factor  in  dispersion,  and  by  the  time 


the  point  of  maximum  ground  level  concentration  is  reached,  the  effect  of 
self-generated  turbulence  on  dilution  is  generally  negligible. 

Combining  Momentum  and  Buoyancy  Contributions  to  Plume  Rise 

Well  defined  methods  are  available  to  predict  the  combined  effects  of 
vertical  jet  momentum  and  buoyancy  on  the  plume  rise  during  its  initial  bent- 
over  phase  when  it  is  influenced  entirely  by  self  generated  turbulence.  How- 
ever, once  the  plume  interacts  with  atmospheric  turbulence  and  levels  off  at  a 
final  rise  height,  there  are  no  methods  for  predicting  this  final  rise  which 
combine  both  the  effects  of  momentum  and  buoyancy.  For  stack  design,  the  usual 
approach  has  been  to  neglect  the  source  momentum  and  compute  only  the  final 
buoyancy  rise.  However,  for  the  high  velocity  release  from  a pipeline  rupture 
or  well  blow-out,  jet  momentum  may  contribute  a significant  fraction  of  the 
total  final  rise.  To  account  for  this,  theoretical  estimates  for  final  rise 
developed  by  Briggs  (1975)  are  extended  to  account  for  both  momentum  and  buoy- 
ancy effects.  The  theory  which  results  is  complicated,  and  a simple  approxima- 
tion is  suggested  in  which  the  separate  momentum  jet  and  buoyant  plume  final 
rises  are  added  using  a sum-of-cubes  law.  This  simple  approximation  is  shown  t( 
yield  accurate  final  rise  predictions  within  ^ 10%  of  the  exact  theoretical 
val ue. 
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CHAPTER  1 


INTRODUCTION 


Plumes  from  high  pressure  pipeline  ruptures  or  well  blowouts  are  consider- 
ably different  from  the  stack  gas  plumes  to  which  most  dispersion  models  are 
tuned.  Some  important  differences  between  analyzing  stack  gas  plumes  and 
dispersion  from  pipeline  ruptures  are: 

• The  release  rate  from  a pressure  vessel  or  pipeline  will  decrease 
with  time  as  the  pressure  in  the  system  drops.  This  will  produce 
a transient  plume  whose  dispersion  characteri sti cs  will  lie  some- 
where between  that  of  an  instantaneous  puff  release,  and  a steady 
plume. 

• The  high  internal  pipeline  pressure  will  cause  the  flow  to  choke 
at  sonic  velocity  at  the  pipe  exit,  and  to  suddenly  expand  and 
accelerate  outside  the  exit  as  it  comes  into  equilibrium  with 
atmospheric  pressure.  Both  the  release  rate  and  outside  expan- 
sion require  the  use  of  compressible  gas  dynamics  to  estimate 
mass  and  momentum  fluxes. 

• The  high  velocity  jet  from  a rupture  will  undergo  considerable 
dilution  due  to  internally  generated  turbulence  before  it  is  bent 
over  by  the  wind  and  dispersed  by  atmospheric  turbulence. 

• Most  accidental  ruptures  occur  at  ground  level,  and  the  plume 
from  a high  velocity  ground  level  source  behaves  differently  than 
a buoyant  plume  from  a high  stack.  Some  method  must  be  found  to 
include  both  the  effects  of  jet  momentum  and  buoyancy  on  the 
trajectory  and  final  rise  of  a high  velocity  pipe  rupture 

rel ease. 

• Concentration  fluctuations  are  present  in  all  plumes  dispersed  by 
turbulence.  However,  it  is  of  crucial  importance  to  accurately 
estimate  these  fluctuations  for  the  highly  toxic  gases  that  may 
be  released  from  pipeline  ruptures  and  well  blowouts  to  determine 
the  risk  of  exposure  to  high  peak  levels.  Because  these  fluctua- 
tions depend  strongly  on  the  source  size  and  conditions  at  the 
point  of  release,  as  well  as  the  intensity  and  scale  of  the 
atmospheric  turbulence  which  disperses  the  plume,  accurate 
estimates  are  much  more  difficult  to  obtain  than  for  mean 
concentrations. 

To  account  for  these  differences  a study  was  carried  out  which  resulted  in  the 
report  "The  Release  and  Dispersion  of  Gas  From  Pipeline  Ruptures"  (January 
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1979).  At  that  time  there  was  little  theory  and  almost  no  experimental  data 
on  which  to  base  a release  rate  and  dispersion  model  for  pipeline  ruptures. 
The  report  recognized  this  and  took  the  point  of  view  that  "it  is  better  to 
make  an  estimate,  which  may  be  at  best  an  educated  guess,  than  to  ignore  the 
effect  of  a parameter  entirely."  The  fundamental  assumptions  on  which  the 
model  was  based  have  remained  supri singly  durable  in  spite  of  the  large 
measure  of  educated  (and  sometimes  desperate)  guesswork  that  went  into  its 
formulation.  A series  of  subsequent  reports;  "Expansion  and  Plume  Rise  of 
Gas  Jets  from  High  Pressure  Pipeline  Ruptures"  (April  1981),  "Predicting  Risk 
of  Exposure  to  Peak  Concentrations  in  Fluctuating  Plumes"  (December  1982)  and 
"Exposure  Time  Effects  on  Concentration  Fluctuations  in  Plumes"  (January 
1985)  have  modified  and  refined  the  theory  and  compared  it  with  both  full 
scale  and  wind  tunnel  data. 

This  report  is  the  fifth  in  this  series.  Its  purpose  is  to  develop 
further  the  analytical  methods  in  the  existing  theoretical  model.  While  its 
general  focus  is  on  the  dispersion  of  toxic  gas  from  pipeline  ruptures,  many 
of  the  theoretical  relations  that  are  developed  can  be  applied  equally  well 
to  other  atmospheric  dispersion  situations.  An  attempt  has  been  made  to 
review  the  existing  results  before  presenting  new  theories  or  recommending 
changes.  Hopefully,  this  will  allow  the  report  to  make  sense  without  the 
need  for  the  reader  to  constantly  refer  to  the  previous  reports  in  the 
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CHAPTER  2 


DISPERSION  FROM  TRANSIENT  RELEASES 

The  release  of  toxic  or  flammable  gases  is  usually  a short  term  event 
that  occurs  over  a time  interval  of  a few  seconds  to  a few  hours.  The 
situation  is  further  complicated  by  release  rates  which  vary  with  time,  for 
example  the  case  of  an  exponentially  decreasing  release  from  the  rupture  of  a 
high  pressure  vessel  or  pipeline.  In  this  chapter  we  will  develop  methods 
for  dealing  with  transient  releases  within  the  traditional  framework  of 
atmospheric  dispersion  of  continuous  steady  plumes  and  instantaneously 
released  puffs. 

Limiting  Cases:  Puffs  and  Steady  Plumes 

A very  short  duration  release  can  be  modelled  as  an  instantaneous  puff, 
while  at  the  other  extreme,  a release  of  long  duration  may  be  thought  of  as  a 
steady  plume,  except  near  its  "nose"  and  "tail"  where  it  starts  and  stops. 
There  are  two  fundamental  differences  between  a puff  and  a plume,  one  related 
to  its  geometry  and  the  other  to  the  way  in  which  its  spread  is  measured.  A 
puff  spreads  in  all  three  dimensions,  x,  y and  z as  it  is  carried  downwind, 
while  a plume  is  diluted  only  by  its  crosswind  spreading  in  y and  z.  The 
alongwind  spreading  of  the  series  of  puffs  that  make  up  a plume  is  ineffect- 
ive in  causing  dilution  because  it  is  material  from  the  puffs  immediately 
upwind  and  downwind  that  is  being  entrained,  rather  than  fresh  air. 

The  ensemble  averaged  dispersion  of  a puff  is  measured  relative  to 
the  mass  center  of  each  of  the  puffs  in  the  ensemble.  Because  the  spread  is 
relative  to  the  mass  centers  of  the  puffs  it  does  not  include  dispersion 
caused  by  meandering  of  the  mass  centers  themselves.  A receptor  fixed  on  the 
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ground  will  see  dispersion  caused  by  both  and  as  it  is  exposed  to  the 
passage  of  an  ensemble  of  puffs.  The  total  crosswind  spread  will  be 
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a.'  + a 

iy  > 
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my 


(2-1) 


where 


a. 


= crosswind  "instantaneous"  spread  of  puffs  to  their  mass  centers 


a = crosswind  "meandering"  spread  of  the  mass  centers  of  the  ensemble 


Because  a continuous  plume  is  just  a series  of  superposed  puffs  the  same 
definitions  apply,  with  the  instantaneous  spread  found  by  taking  an  ensemble 


instantaneous  concentration  profile.  The  total  dispersion  is  found  either  by 
taking  a long  time  exposure,  or  by  measuring  the  meandering  of  the  centroids 
of  each  ensemble  of  snapshots  and  using  (2-1). 

The  important  distinction  between  the  spread  of  puffs  and  plumes  is 
simply  one  of  terminology.  The  term  "spread"  used  in  the  literature  is 
almost  always  used  to  refer  to  the  instantaneous  "relative"  dispersion  of  a 
puff  and  to  the  total  dispersion  of  a continuous  plume.  Unless  this  differ- 
ence is  understood  a great  deal  of  confusion  can  result  in  interpreting  puff 
and  plume  spreads  reported  in  the  literature.  From  the  standpoint  of  risk 
analysis,  we  are  interested  in  a fixed  receptor  on  the  ground,  and  the 
dispersion  of  an  ensemble  of  puffs  is  the  same  as  that  of  a plume,  and 
includes  both  the  instantaneous  and  meandering  spreads. 

The  geometrical  differences  between  puff  and  plume  dispersion  appear 
explicitly  in  the  equations  for  maximum  concentration  at  the  center  of  an 
ensemble  of  puffs,  each  of  which  has  a total  mass  release  of  m^,  kg. 


of  puffs 
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(271)  a a a 
X y z 


For  the  same  mass  released  slowly  over  a long  time  duration 
may  be  modelled  as  a continuous  plume,  with 

c Q 

plume  2it  U o a 
c y z 


(2-2) 

t^  the  dispersion 
(2-3) 


where  Q = is  the  source  strength.  Taking  the  ratio  of  (2-3)  to  (2-4) 

and  defining  a release  length  x = t U , where  U is  the  mean  convection 
velocity  of  plume  material  through  a plume  crosswind  section,  we  obtain 

(!^)  (2-4) 

C ..  X ^ 

puff  r 
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The  ratio  is  unity  when  a = (2ti)  * x which  is  a - 0.4  x . For  short 

duration  releases,  where  x^  is  small,  the  plume  model  predicts  higher 
concentrations  than  the  puff  model  because  the  plume  neglects  the  alongwind 
dispersion.  For  long  duration  releases  the  puff  model  yields  higher 
concentrations  because  all  the  mass  m^  is  released  instantaneously 
rather  than  spread  out  over  a distance  x^. 

Transition  From  Puff-like  to  Plume-like  Dispersion 

From  the  preceding  discussion  it  is  clear  that  at  the  two  limits  of 
short  puff-like  bursts  and  long  plume-like  releases,  the  correct  model,  (2-2) 
or  (2-3)  is  the  one  which  predicts  the  lowest  concentration.  That  is,  as  a 
model  becomes  i nappropri a te  it  will  tend  to  overestimate  the  concentration. 
This  observation  provides  a simple  way  of  choosing  whether  a release  is 
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puff-like  or  plume-like;  simply  evaluate  the  concentration  using  both  models 
and  choose  the  one  that  yields  the  smallest  concentration. 

Another  approach  to  defining  a transition  between  puff  and  plume-like 
behavior  is  to  examine  how  the  upwind  "nose"  and  downwind  "tail"  of  a 
transient  release  modify  the  continuous  plume  between  them.  This  situation 
is  shown  schematically  in  Figure  2-1  where  the  release  is  viewed  from  above 
after  having  travelled  some  distance  downwind.  The  concentration  distribu- 
tion has  a nose  and  tail  that  are  influenced  by  and  disperse  like  puffs. 
Between  them  is  a region  which  behaves  like  a continuous  plume,  uninfluenced 
by  the  additional  dilution  of  the  nose  and  tail.  As  the  release  moves  down- 
wind, increases  while  remains  constant,  and  the  nose  and  tail  puffs 

gradually  erode  the  plume-like  core.  For  simplicity,  we  will  assume  that  the 

-0  5 

plume-like  core  region  disappears  when  = {2%)  * x^,  which  is  the  point 

where  the  plume  and  puff  models  in  (2-4)  predict  the  same  concentration. 

From  Figure  2-1  we  see  that  the  nose  and  tail  merge  when  R = 0.5  x^,  so  that 

0 5 

the  effective  "radius"  of  these  puffs  must  be  R = (tc/2)  * a^,  which  seems 

reasonabl e. 

The  transition  to  complete  puff -like  behaviour  requires  the  nose  and 

tail  zones  to  completely  merge.  Based  on  the  experiments  of  Fackrell  and 

Robins  (1981)  for  ground  level  sources,  the  transition  to  puff-like  behavior 

is  complete  when  the  nose  (or  tail)  puff  radius  R has  increased  to  about 

0 5 

three  times  the  critical  value  {n/2)  * where  the  nose  and  tail  puffs  begin 
to  merge.  Using  this  observation,  the  criteria  for  release  behavior  are: 


a < 0.4  X 
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pi  ume-1 i ke 


a^>  1.2  X 
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(2-5) 


Distributed  mass  flux 
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Fig.  2 


-1  Along-wind  diffusion  at  the  nose  and  tail  of  a transient  suddenly 
started  and  stopped  plume. 
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A release  of  fixed  length  will  behave  initially  as  a plume,  until 

along-wind  diffusion  of  the  nose  and  tail  cause  a transition  to  puff-like 
behaviour.  Wilson  (1981a.)  gives  a more  rigorous  treatment  of  diffusion  in 
the  nose  and  tail  of  a plume,  and  develops  a plume  model  that  includes  an 
along-wind  component. 

The  criteria  in  (2-5)  are  in  good  agreement  with  the  wind  tunnel  data  of 
Fackrell  and  Robins  for  diffusion  of  a large  ensemble  of  100  to  300  releases 
from  a ground  level  source.  These  experiments  also  confirm  that  during  the 
transition  from  plume-like  to  puff-like  behaviour,  both  puff  and  plume  models 
overestimate  the  concentration.  The  puff  model  overestimates  concentration 
because  it  concentrates  the  release  at  a single  instant,  and  the  plume  model 
also  overestimates,  because  it  neglects  along-wind  dispersion.  Thus, 
selecting  the  smaller  of  the  two  predictions  will  always  be  conservative. 

Characteristic  Release  Duration 

The  previous  discussion  considered  the  simplest  type  of  transient 
release;  a suddenly  started  and  stopped  steady  release  rate.  Most  practical 
situations,  particul arily  accidental  releases  from  pipe  and'  pressure  vessel 
ruptures,  have  a time-varying  release  rate  m^.  The  choice  of  an  appropriate 
release  duration  is  not  obvious,  particul ari ly  for  releases  which  start  and 
stop  for  short  periods  during  the  release.  Robins  (1984)  suggests  that  the 
best  choice  is  the  mass-weighted  time  centroid  defined  as 

t.  = — / m t dt  (2-6) 

0 ® 

where  m^  is  the  total  mass  released 
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m = J m dt  (2-7) 

r ^ e ^ 

0 


For  an  exponentially  decaying  release  rate  from  a pipe  rupture 


it  is  easy  to  show  that  (2-8)  yields 


so  that  the  characteristi c release  duration  is  equal  to  the  decay  time 
constant. 


Estimating  Along-Wind  Spread 

Predicting  transition  from  plume-like  to  puff-like  behaviour  requires  an 
estimate  for  the  along-wind  spread  This  component  of  plume  spread  is 

strongly  influenced  by  vertical  wind  shear  dU/dz,  which  enhances  the  rate  of 
spread.  As  the  plume  spreads  vertically  wind  shear  causes  an  along-wind 
separation  of  the  plume  layers,  which  are  then  smeared  together  by  further 
vertical  spreading.  Wilson  (1981a)  shows  that  for  a ground  level  source  in  a 
power  law  velocity  profile  with  U that  the  shear  induced  along-wind 

spread  is 


^shear 


0.6  n (2.9)"  X 


SO  that  the  total  along-wind  spread  is 

2 2 

o = a 

^ ^turbulence 


2 

a 

^shear 


(2-10) 


(2-11) 


The  strong  x dependence  of  the  shear  component  usually  makes  it  the  dominant 


term. 
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It  is  interesting  to  note  that  for  o^<^x  the  wind  shear  component  in 

3 /2 

(2-10)  predicts  that  x ' . This  is  the  same  exponent  predicted  for 
relative  dispersion  of  puffs  in  the  inertial  subrange  of  turbulence,  and  some 
of  the  data  reported  by  Gifford  (1977)  which  shows  this  accelerated  spreading 
may  in  fact  be  dominated  by  wind  shear  rather  than  turbulence. 


Estimating  Instantaneous  Plume  or  Puff  Spread 

In  the  previous  discussion  the  total  spread  in  (2-1)  was  shown  to  be 
the  appropriate  parameter  to  describe  the  dispersion  of  an  ensemble  of  puffs 
or  instantaneous  plumes  relative  to  a specific  fixed  receptor  on  the  ground. 
However,  only  the  instantaneous  spread  represents  real  dilution  of  source 
material  by  entrainment  of  ambient  air.  The  meandering  component  is 
generated  by  large  scale  atmospheric  turbulence  flapping  the  instantaneous 
plume  (or  a series  of  puffs)  back  and  forth  over  the  receptor.  This  random 
motion  of  the  plume  or  puff  mass  center  causes  an  apparent  dispersion  to  a 
fixed  observer,  as  the  plume  material  meanders  away  to  contaminate  some  other 
receptor  location. 

In  a situation  where  there  is  no  specific  receptor  location,  for  example 
in  estimating  the  general  toxicity  or  flammability  hazard  at  a given  radius 
from  a source,  it  is  more  useful  to  use  the  instantaneous  plume  spreads 
and  o.^  as  measures  of  dilution.  Because  atmospheric  diffusion  models  used 
in  design  and  regulation  of  emissions  predict  only  the  total  plume  spreads, 
an  estimate  of  instantaneous  spread  in  terms  of  the  total  spread  is  the  most 
useful  form  in  which  to  express  A detailed  discussion  of  functional 

relationships  for  total  and  instantaneous  spread  is  given  in  Appendix  A where 
it  is  shown  that  total  spread  responds  to  the  turbulence  spectrum  with  high 
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frequency  components  filtered  out,  and  instantaneous  spread  to  a band-pass 
filter  that  removes  both  high  and  low  frequency  components.  Roth  and 
share  a dependence  on  the  crosswind  turbulence  intensity  v'/U^,  the  Eulerian 
turbulence  scale  and  its  Lagrangian  counterpart,  L^. 

In  terms  of  the  turbulence  intensity  and  scale,  the  power  law  equation 
for  crosswind  spread  is  given  by  (A2) 


2 2ay  2 
D X + a 

y s 


-,0.5 


with  the  constant  given  by  (A8) 


D = (2A,  A 
v ^ 1 


K] 


(2-12) 


(2-13) 


with  a value  of  A^  = 0.6  recommended  by  Pasquill  (1974),  pp.  89-92. 

In  this  equation  the  effect  of  source  size  is  added  in  the  form  of  a 
random  uncorrelated  variance  to  the  turbulent  spread.  This  random 
contribution  of  source  size  is  a consequence  of  a fundamental  assumption  in 
the  statistical  theory  of  dispersion,  which  simply  tags  fluid  particles  as 
they  pass  the  source  location.  These  tagged  "source"  particles  immediately 
take  on  the  turbulence  properties  of  the  diffusing  flow  and  account  for 
source  size  only  as  a random  variable  on  initial  particle  position.  This 
leads  to  source  size  appearing  as  an  additive  variance  in  (2-12). 

Most  real  sources  are  formed  by  a fluid  jet,  which  remains  coherent 
after  release.  Some  distance  after  release  is  required  before  the  turbulence 
velocities  and  scales  of  the  diffusing  flow  become  dominant  and  decorrelate 
the  internal  flow  structure  of  the  source  jet.  For  real  sources  a linear 
addition  of  source  size  may  be  a better  approximation,  and  (2-12)  takes  on 
the  form 
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dw 

a = D X + a 

y y s 


(2-14) 


If  the  instantaneous  plume  spread  is  dominated  by  turbulence  in  the  inertial 
subrange,  it  is  predicted  by  (A28)  in  Appendix  A 
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which  for  = 0 reduces  to  the  expected  relationship,  x . This 

equation  is  valid  only  for  spreading  close  to  the  source  where  the  plume  or 
puff  width  is  small  compared  to  the  turbulence  scale  In  these 

ci rcumstances  the  total  spread  grows  linearly  with  x,  so  that  a^  = 1.0. 
Using  this  value  in  (2-13)  and  (2-14) 


a = ( 1-  ) X + a 

y u s 

c 


(2-16) 


near  the  source.  Using  (2-16)  in  (2-15)  the  ratio  of  instananeous  to  total 
spread  is  given  by 
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(2-17) 
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which  is  valid  for  a.  <<  A . For  the  full  turbulence  spectrum,  the  result 
iy  V K > 

is  more  complicated.  Using  (2-16)  in  (A40) 
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(2-18) 
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which  should  be  valid  as  long  as  < A^. 
from  (A-39) 


The  source  size  factor  is, 


$ 

s 


(2-19) 


Experimental  Verification  of  Instantaneous  Spread  Predictions 

Although  there  is  a great  deal  of  experimental  data  on  total  and 
instantaneous  plume  spreads,  only  the  data  of  Nappo  (1984)  presents  measure- 
ments of  both  quantities  in  the  same  flow,  along  with  the  turbulence 
statistics  needed  to  test  the  predictions  of  (2-17)  and  (2-18).  Nappo  used 
photography  to  record  an  ensemble  of  16  instantaneous  snapshots  of  smoke 
plume  dispersing  in  grid  turbulence  in  a wind  tunnel.  These  16  independent 
realizations  used  to  measure  instantaneous  spreads  were  taken  at  15  second 
intervals  during  which  a 4 minute  time  exposure  was  recorded  to  measure  total 
spread.  Nappo' s objective  was  to  deduce  turbulence  intensity  and  scale  from 
the  observed  smoke  plume  spreads.  Here  we  will  use  the  measured  spreads  and 
turbulence  statistics  from  hot  wire  anemometer  measurements  to  test  the 
equations  for  estimating  instantaneous  spread. 

The  turbulence  grid  used  by  Nappo  was  a square  bar,  square  mesh  biplane 
grid  with  a mesh  size  of  10  cm  and  a bar  size  of  2 cm.  Square  plates 
30  X 30  cm  were  used  to  block  9 adjacent  holes  out  of  every  36  to  create  high 
intensity  large-scale  turbulence.  A circular  tube  with  an  inside  diameter  of 
0.873  cm  injected  smoke  isoki neti cal ly  into  the  flow  400  cm  downstream  from 
the  grid.  Measurements  of  vertical  spreads  and  were  made  by  Nappo 
from  photographs  at  8 locations  from  x = 25  cm  to  x = 350  cm. 
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There  are  were  two  major  obstacles  in  comparing  Nappo's  measurements 
with  the  theory.  The  first  is  the  natural  decay  of  grid  turbulence  with  down- 
stream distance.  The  x direction  intensity  and  scale  varied  from  u'/U  = 0.13 
and  - 21  cm  at  the  source  to  0.07  and  33  cm  at  x = 400  cm.  In  contrast, 

the  theoretical  predictions  assume  homogeneous  turbulence  with  constant  in- 
tensity and  scale.  The  second  obstacle  is  that  only  x direction  turbulence 
statistics  were  measured,  using  a single  sensor  hot  film  anemometer.  This 
required  the  cross-stream  statistics  to  be  estimated  from  along  wind  values. 
Grid  turbulence  is  generally  anisotropic  in  intensity  but  has  integral  scale 

relations  that  are  close  to  isotropic.  For  grid  turbulence  we  estimated 
2 2 2 

that  u‘  ==  2w‘  - 2v'  , and  that  along-wind  and  crosswind  scales  were 

isotropically  related,  with  = 0.5  A^.  With  the  scales  estimated  by 
fitting  exponentials  to  Nappo's  autocorrel a tions,  the  average  intensity  and 
scale  over  the  350  cm  plume  length  was  found  to  be  v'/U  - 0.067  and  A^  - 16 
cm. 

Nappo's  measurements  for  total  plume  spread  were  used  to  obtain  a 
refined  estimate  of  turbulence  intensity  and  effective  source  size.  The  data 
are  shown  in  Figure  2-2  where  they  are  in  good  agreement  with  (2-16)  with 
v'/U^  = 0.062  = 0.70  cm,  and  = 1.0  for  grid  turbulence  as  measured  by 

Snyder  and  Lumley  (1971).  Because  intensity  and  scale  are  inversely  related 
in  grid  turbulence,  this  lower  effective  intensity  implies  a larger  effective 
scale  see  (A6),  and  A^  = 17  cm  was  taken  as  the  refined  estimate. 

Figure  2-2  shows  the  predictions  of  (2-17)  and  (2-18)  using  these  values 
for  intensity,  scale  and  source  size.  The  strong  effect  of  source  size  on 
the  instantaneous  spread  is  evident,  as  well  the  persistent  difference 
between  instantaneous  and  total  spreads.  In  general  the  experimental  results 
are  in  good  agreement  with  the  full  spectrum  predictions. 
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Figure  2-2  Instantaneous  a-jy  (triangles)  and  time-averaged  Cy  (circles)  spread  of  a plume 
from  a d = 0.7  cm  source  in  grid  turbulence.  Data  from  Nappo  (1984).  Theory 
from  (2-15),  (2-16)  and  (2-18),  (2-19)  with  Ai  = 1.0,  v'/Uc  = 0.062  Av  = 17  cm 


Other  theoretical  estimates  are  available  for  the  instantaneous  to  total 
plume  spread  ratio.  Hanna  (1984)  used  Gifford's  (1982)  Langevin  equation 
formulation  to  solve  for  ratio  of  instantaneous  to  total  spread. 


where 


0.5  S^(l  - exp(-  T')) 


(2-20) 


2A! 


+ r - (1  - exp(-  T')) 


and  = 0.6  is  the  Eul  eri an-Lagrangi an  scale  ratio  constant  in  (A6).  For 
most  practical  situations  where  < 0.1,  Hanna  shows  that  - 1.0  is  a 

good  approximation. 

Figure  (2-3)  shows  a comparison  of  Hanna's  equation  and  our  theory  from 
(2-17)  and  (2-18)  with  Mappo's  experiments.  As  expected,  the  approxima tion 
of  the  spectrum  as  an  inertial  subrange  (-  5/3  spectrum)  applies  only  very 
close  to  the  source  where  <<  A^.  Hanna's  theory  predicts  much  smaller 
values  of  the  instantaneous  to  total  spread  ratio  near  the  source.  In  part 
this  is  due  to  the  use  of  an  effective  source  size  derived  from  a linear 
source  correction  in  (2-14).  If  an  additive  variance  is  assumed,  using 
(2-12),  we  find  the  total  spread  is  well  represented  by  v'/U^  = 0.067  and 
= 1.4  cm.  With  these  values  and  = 16  cm  Hanna's  model  in  (2-20) 
performs  better,  predicting  a minimum  variance  ratio  of  0.20  at  x = 60  cm, 
close  to  the  observed  value. 
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Figure  2-3  Ratio  of  instantaneous  to  total,  time  averaged,  plume  spread  variances  in  grid 
turbulence  experiments  of  Nappo  (1984).  Theory  from  (2-15),  (2-16)  and 
(2-18),  (2-19).  Hanna's  theory  from  (2-20);  a = 0.7  cm.  A,  = 1.0. 
v'/Uc  = 0.062,  Av  = 17  cm  ^ ^ 
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We  conclude  that  either  model  gives  an  adequate  representation  of  the 
ratio  of  instantaneous  to  total  spread.  Because  Hanna's  model  always  reduces 
to  the  correct  asymptotic  value  of  unity  near  to  and  far  from  the  source,  it 
is  the  preferred  formulation.  On  the  other  hand,  the  model  in  (2-17)  gives  a 
better  physical  insight  into  the  influence  of  the  turbulence  spectral 
bandpass  on  instantaneous  spread. 

Summary 

In  this  chapter  a set  of  criteria  have  been  developed  to  mark  the 
transition  of  a release  from  its  initial  plume-like  behaviour  to  puff-like 
dispersion  farther  downwind.  Because  a plume  dispersion  model  using  will 
always  overestimate  the  concen tration  through  its  neglect  of  along-wind 
diffusion,  the  use  of  plume  models  for  transient  releases  will  give  conserva- 
tive estimates  of  exposure  hazards. 

Dispersion  relative  to  the  instantaneous  center  of  mass  of  a puff  or 
plume  concentration  profile  is  shown  to  be  a useful  way  of  character! zi ng  the 
actual  dilution  of  source  material.  Relationships  between  this  instantaneous 
spread,  and  the  total  spread  which  contains  the  meandering  component  of 
dispersion  are  developed.  By  relating  instantaneous  to  total  plume  spread 
the  readily  available  total  spread  relations  in  each  atmospheric  stability 
class  can  be  directly  applied  to  hazard  assessment. 


( 


CHAPTER  3 


SOURCE  SIZE  EFFECTS  ON  CONCENTRATION  FLUCTUATIONS 

The  most  remarkable  charac teri Stic  of  concentration  fluctuation 
intensity  is  its  sensitivity  to  the  geometry  and  release  conditions  at  the 
source.  To  an  engineer  whose  experience  dictates  that  the  influence  of  source 
geometry  on  a jet  is  limited  to  distances  less  than  SO  source  diameters,  it 
seems  incredible  that  concen trati on  fluctuations  will  be  influenced  by  source 
size  at  distances  several  kilometers  downwind.  However,  for  small  elevated 
sources  dispersing  in  the  large  scale  atmospheric  turbulence  above  the  sur- 
face layer,  this  strong  sensitivity  is  certainly  present,  and  increases  the 
difficulty  in  predicting  peak  concentration  levels. 

The  physical  mechanism  behind  this  source  size  sensitivity  lies  in  the 
difference  between  instantaneous  plume  spread  and  meandering.  The  source 
size  acts  as  a filter  which  determines  the  amount  of  small  scale  atmospheric 
turbulence  that  will  be  effective  in  diluting  the  plume  by  entrainment.  The 
larger  scales  do  not  contribute  to  air  entrainment  but  instead  cause  the 
narrow  instantaneous  plume  to  meander  back  and  forth.  This  generates  a high 
level  of  i n termi ttency  and  produces  a high  concentration  fluctuation  intens- 
ity at  a fixed  receptor. 

The  instantaneous  plume  spread  from  a small  source  increases  very  slowly 
in  large  scale  turbulence,  because  it  filters  out  most  of  the  turbulent 
velocity  fluctuations.  Because  the  width  of  an  instantaneous  plume  from  a 
small  source  remains  small,  it  is  subject  to  meandering  for  much  longer 
distances  than  the  plume  from  a large  source.  The  details  of  this  turbulence 
filtering  are  discussed  in  Appendix  A. 
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When  periods  of  zero  concentration  are  neglected,  and  conditional  mean 
and  fluctuation  statistics  are  measured,  the  effect  of  source  size  can  gener- 
ally be  neglected,  as  shown  by  Wilson  (1982)  and  Wilson,  Robins  and  Fackrell 
(1985).  The  reason  for  this  change  from  strong  source  size  sensitivity  for 
total  variance  to  almost  complete  insensitivity  for  conditional  variance  is 
that  meandering  affects  the  plume  by  producing  intermittent  periods  of  zero 
concentra tion.  When  these  are  removed  in  computing  the  conditional  statis- 
tics, the  remaining  small  scale  mixing  process  is  much  less  dependent  on 
source  size. 


Along-Wind  Variation  of  Fluctuation  Intensity 

Even  in  homogeneous  turbulence  far  from  the  influence  of  ground  surface 

effects,  concentration  fluctuation  intensity  will  vary  with  downwind 

distance.  This  variation  is  due  to  the  balance  between  advection,  produc- 

2 

tion,  diffusion  and  dissipation  of  the  total  variance  c as  well  as  the 
decrease  in  the  mean  concentration  “c.  By  using  the  experimentally  observed 
balance  between  advection  and  dissipation  of  variance,  and  a simple  model  for 
the  dissipation  time  scale,  Wilson,  Fackrell  and  Robins  (1982)  derived  a 
relation  for  the  total  fluctuation  intensity  on  the  plume  axis. 


0.34X 

.2  ,2 

\ + X 

0 


(3-1) 


where  i is  the  total  fluctuation  intensity  on  the  plume  axis  in  the  absence 
of  wind  shear  or  ground  surface  effects.  By  expressing  the  along-wind 
distance  \ = in  plume  coordinates,  the  theory  attempts  to  provide  a 

universal  formulation  that  may  be  valid  for  a range  of  atmospheric  stability 
and  surface  roughness  conditions. 
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The  virtual  origin  \ accounts  for  the  effect  of  source  size  on  the 
0 

total  variance,  and  is  determined  by  the  ratio  of  source  size,  d,  to  the 

average  integral  scale  of  turbulence  A = + A^)/2.  Wilson,  Fackrell  and 

Robins  (1982)  included  this  scale  by  using  the  ratio  of  source  height  h to 

diameter  d in  their  equation  for  the  virtual  origin  For  the  experimental 

data  set  of  Fackrell  and  Robins  (1982a)  used  to  determine  constants  in  the 
theory,  Robins  (1979)  reported  values  of  these  turbulence  scales.  This  data 
may  be  approximately  fit  by  A ==  0.6z,  for  z < 0.25  H,  where  H is  the 
thickness  of  the  atmospheric  shear  layer.  Using  this,  the  empirical  equation 
for  the  virtual  origin  is 


\ 

° 1 + 0.055  L 

d 


(3-2) 


To  account  for  inhomogeneous  shear  flows,  the  scale  A is  evaluated  at  a 

reference  height.  In  (3-1)  and  (3-2)  the  constants  proposed  by  Wilson, 

Fackrell  and  Robins  (1982)  are  adjusted  by  the  factor  1.177  to  account  for 

the  use  of  plume  standard  deviations  /(a  a 1 rather  than  half  widths  /f6  6 ] 

as  the  distance  parameter.  Also,  their  original  formulation  expressed  the 

along-wind  variation  in  terms  of  a variance  source  strength  ratio  q/0  and  a 

separation  factor  p (or  <}))  for  the  pair  of  variance  sources.  Roth  these 

factors  are  combined  in  the  constant  0.34  in  (3-1). 

Sykes,  Lewellen  and  Parker  (1984)  propose  a model  for  the  dissipation 

time  scale  that  accounts  di rectly  for  the  effect  of  source  size  on  total 

variance.  In  the  appendix  to  their  paper  they  derive  a simple  model  for  i 

2/3  1/3  “1 

that  predicts  the  distance  X to  maximum  intensity  as  X d A H . 

This  result  is  a consequence  of  limiting  the  turbulence  to  the  inertial  sub- 
range of  the  turbulence  spectrum  where  the  kinetic  energy  of  turbulence 
decays  as  the  -5/3  power  of  frequency. 
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Sawford  (1983a)  (1983b)  using  random-walk  Lagrangian  particle  statistics 

shows  that  it  is  the  Lagrangian  rather  than  the  Eulerian  turbulence  scale  that 

is  important.  Using  the  Eulerian  to  Lagrangian  scale  relation  in  equation 

(A6)  from  Appendix  A,  his  results  predict  the  distance  to  maximum  vari- 

2/3  1/3  -1 

ance  should  be  proportional  to  d ^ A ^ where  u^  is  the  RMS  turbulent 

velocity  scale.  It  is  shown  in  Appendix  R that  the  turbulence  velocity  scale 

u /U  in  X disappears  when  the  plume  spread  coordinate  \ is  used, 
s c m3  X m3  X 

In  contrast,  the  empirical  equation  (3-1)  of  Wilson,  Fackrell  and  Robins 

(1982)  has  a maximum  intensity  at  X = X . For  small  source  sizes  where  A/d 

is  large,  (3-2)  yields  X = 2.2  d a"\  so  that  X d A~^.  This  empirical 

^ 0 max  ^ 

fit  to  the  wind  tunnel  data  is  limited  to  the  range  4 < A/d  < 40.  The  differ- 

2/3 

ence  between  the  theoretical  prediction  of  X « d and  the  empirical 

^ max 

equation  wi  th  X « d is  significant. 

^ max  ^ 

The  greatest  weakness  of  the  virtual  origin  correction  in  (3-2)  is  that 

it  is  only  an  empirical  correlation,  tuned  to  one  specific  neutrally  stable 

-4 

wind  tunnel  shear  flow,  with  a surface  roughness  of  z^/H  = 4 x 10  , (about 
z^  = 25  cm  in  the  atmosphere).  In  the  next  section  we  will  attempt  to  reduce 
this  deficiency  by  deriving  relations  for  source  size  directly  from  turbulence 
parameters. 

Source  Size  Effects  Derived  from  Turbulence 

The  major  contribution  of  Sykes,  Lewellen  and  Parker  (1984)  to  develop- 
ing an  improved  theory  for  the  along-wind  development  of  concentration 

2 

fluctuations  was  their  suggestion  that  the  dissipation  time  scale  for  c’ 
should  be  governed  by  variables  that  influence  the  instantaneous  plume  spread, 
rather  than  the  total  plume  spread  which  includes  meandering.  Because 
meandering  moves  the  entire  plume  back  and  forth,  it  will  not  contribute  to 


dissipation,  which  is  governed  by  processes  that  cause  actual  dilution  of 
pi  ume  ma teri al . 
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They  developed  a scaling  hypothesis  for  the  portion  u^^.  of  the  total  RMS 
turbulence  velocity  u^  that  is  effective  in  contributing  to  the  spread  of 

the  instantaneous  plume,  and  suggested 

a.  1/3 


This  velocity  scaling  implies  an  inertial  subrange  turbulence  spectrum  with 
kinetic  energy  proportional  to  the  -5/3  power  of  fluctuation  frequency. 

In  Appendix  B these  scaling  principles  are  used  to  develop  a simple 
theory  for  along-wind  variation  of  total  fluctuation  intensity.  In  homoge- 
neous turbulence,  with  constant  integral  scale  and  RMS  velocity  the  theory 
predicts  in  (B22) 


(3-4) 


The  effect  of  source  size  is  contained  in  the  virtual  origin  For 

turbulence  in  the  inertial  subrange,  the  effect  of  source  size  and  turbulence 
scale  on  is  given  by  (B23).  Because  the  source  size  must  be 
proportional  to  source  diameter,  (B23)  can  be  written  as 


X 

s 


0 


(3-5) 


For  homogeneous  turbulence,  where  A is  a constant,  the  maximum  intensity  in 
(3-4)  occurs  at 


X 


X 


(a-1) 


max 


(3-6) 
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at  which  point  the  maximum  intensity  on  the  plume  axis  is 


The  distance  \ to  maximum  concentration  from  (3-5)  and  (3-7)  is 
max  ^ ^ ' 

2/3  1/3  -1 

proportional  to  d ' A ' H which  is  in  good  agreement  with  the  wind  tunnel 
data  reported  by  Fackrell  and  Robins  (1982b)  for  an  elevated  source  in  a weak 
shear  flow.  The  data  of  Fackrell  and  Robins  (1982a)  has  a dependence  on 
source  size  that  indicates  i « d , which  implies  from  (3-5)  and  (3-8) 

that  a = -1.5. 

Accounting  for  Wind  Shear  and  Inhomogeneous  Turbulence 

The  Wilson,  Fackrell  and  Robins  model  in  (3-1)  and  (3-2)  and  the  new 
theory  of  (3-4)  and  (3-5)  were  developed  for  homogeneous  turbulence  and  no 
wind  shear.  Sawford's  (1983b)  Lagrangian  particle  tracking  simulations 
identified  two  different  effects  of  boundary  layer  turbulence  on  concentra- 
tion fluctuations. 

• As  the  plume  moves  downwind  it  will  sample  varying  scales  of  turbu- 
lence. These  varying  scales  will  change  the  meandering  process  of 
the  instantaneous  plume,  and  al ter  the  effect  of  scale  on 

• The  mean  wind  shear  dU/dz  stretches  and  distorts  eddies,  increasing 
their  mixing  rates  and  causing  the  fluctuation  intensity  to  rapidly 
lose  its  dependence  on  source  size. 

Wilson  (1982)  suggests  a simple  correction  to  the  theory  to  simultaneously 
account  for  effects  of  changing  scales  and  of  wind  shear  on  the  plume  from  an 
elevated  source,  emitted  at  z = h in  a boundary  layer.  These  effects  are  ac- 
counted for  by  evaluating  the  scale  A used  in  (3-2)  at  a reference  height  z^ 
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= h exp 


(3-8) 


with  A = 0. 6 in  (3-2)  for  the  neutrally  stable  wind  tunnel  boundary  layer 


becomes  independent  of  source  size. 

Because  this  simple  correction  attempts  to  account  for  wind  shear  and 
varying  turbulence  scale  with  a single  parameter,  it  is  limited  to  shear 
flows  where  the  turbulence  scales  and  wind  shear  are  linked  together.  This 
restriction  may  limit  the  model's  application  to  near-neutral  atmospheric 
stability.  Sawford's  (1983b)  calculations  suggest  that  an  additional  length 
scale  proportional  to  (H  U”^  dU/dz)  may  be  required  to  describe  more  complex 


For  the  inertial  subrange  model  in  (3-4)  and  (3-5)  a single  scale 
adjustment  will  also  be  used  to  account  for  both  turbulence  inhomogeniety  and 
wind  shear  decorrel a tion.  The  wind  tunnel  data  of  Fackrell  and  Robins 
(1982a),  used  to  set  the  model's  constants,  has  a scale  variation  with  height 
of  A = 0.6z.  For  a ground  level  source,  an  appropriate  choice  of  height  at 
which  to  evaluate  the  effective  scale  is  z^^^  = while  for  an 

elevated  source  far  from  the  influence  of  surface  reflection  and  dissipation 
effects  z = h is  a logical  choice.  A simple  interpolation  equation  between 
these  limits  is 


Using  the  wind  tunnel  data,  one  set  of  contants  that  is  physically  realistic 
is:  a = -1.5,  B^  = 0.8,  = 2.5  and  B^  = 0.5.  Using  these,  (3-4)  and  (3-5) 


data  in  the  study.  For  a ground  level  source  A = 0 from  (3-8),  and  (3-2) 


fl ows. 


= [h^  + 


0.5 


(3-9) 


become 
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X = 0.8 
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1/3 
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with  the  reference  scale  evaluated  from  (3-9) 


(3-10) 


(3-11) 


z ^ = h [l 
ref 


0.25( 


2 0.5 


(3-12) 


Predictions  using  these  equations  are  compared  to  the  wind  tunnel  data  in 
Figure  3-1.  Agreement  is  not  as  good  as  the  empirical  fit  shown  in  Figure 
3-2  using  equations  (3-1)  and  (3-2)  with  the  effective  scale  from  (3-8). 
However,  the  assumptions  used  to  derive  (3-10),  (3-11)  and  (3-12)  have  a 
stronger  physical  basis  than  the  purely  empirical  fit,  and  should  be  more 
reliable  when  extrapolated  to  other  atmospheric  stabilities,  surface  rough- 
ness and  source  conditions.  (Note  that  the  predictions  in  Figures  3-1,  3-2, 
include  ground  surface  interaction  effects  discussed  in  the  following 
secti on. ) 

Both  the  empirical  correlation  model  (3-1)  and  (3-2)  and  the  inertial 
subrange  formulation  of  (3-10)  and  (3-11)  predict  a constant  asymptotic  value 
for  the  intensity  from  a ground  level  source,  in  agreement  with  the  measure- 
ments of  Fackrell  and  Robins  (1982a).  However,  this  agreement  is  forced  by 
the  way  in  which  the  turbulence  scale  relations  (3-8)  and  (3-12)  were  chosen 
and  should  not  be  regarded  as  a verification  of  the  theory. 

In  homogeneous  turbulence  where  A and  are  constant  both  models 

predict  a decay  rate  i « far  downwind.  This  result  is  in  striking 
contrast  to  the  Lagrangian  particle  tracking  simulations  of  Sawford  (1983a) 
who  predicts  a constant  asymptotic  value  of  i in  this  situation. 


Total  Fluctuation  Intensity  y,^  Intermittency  at  Plume  Height 
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Fig.  3-1  Inertial  subrange  theory  for  intermittency  and  total  concentration 

fluctuation  intensity  for  an  elevated  source  at  h/H  = 0.192  in  a wind 
tunnel  boundary  layer  of  H = 120  cm.  Data  of  Fackrell  and  Robins 
(1982a)  theory  from  (3-10),  (3-11),  (3-12)  and  (3-13),  (3-14),  (3-15) 
with  a = 0.6,  4)  = 0.706 


Total  Fluctuation  Intensity  Intermittency  at  Plume  Height 
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Fig.  3-2  Semi -empi ri cal  theory  for  intermittency  and  total  concentration 

fluctuation  intensity  for  an  elevated  source  at  h/H  = 0.192  in  a wind 
tunnel  boundary  layer  of  H = 120  cm.  Data  of  Fackrell  and  Robins 
(1982a),  theory  from  (3-1),  (3-2)  and  (3-13),  (3-14)  and  (3-15)  with 
a = 0.6,  (})^  = 0.  706. 
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Unfortunately  there  are  no  experiments  in  homogeneous  turbulence  with  which 
to  resolve  this  difference.  Grid  turbulence,  which  is  often  used  as  an 
"approximately"  homogeneous  flow  actually  has  a length  scale  A which 
increases  with  x,  and  behaves  more  like  A <==  X than  A = constant. 


Surface  Effects  on  Variance  Profiles 

Wilson  (1982)  approximates  the  spatial  variation  of  both  the  mean  "c  and 
2 

variance  c'  profiles  by  Gaussian  profiles.  The  effect  of  wind  shear  is 
modelled  by  using  a single  variance  source  displaced  above  the  actual  source 
height  by  an  amount  where  4>^  = 0.706  was  determined  experimentally. 

The  original  theory  in  Wilson,  Fackrell  and  Robins  (1982)  and  Wilson  (1982) 
used  a pair  of  variance  sources,  which  were  located  on  either  side  of  the 
point  of  release.  However,  Wilson  and  Simms  (1985)  later  showed  a single 
source  gave  more  realistic  predictions  in  the  tails  of  profiles. 

Near  ground  level,  the  presence  of  small  scale  turbulence  and  wind  shear 
increase  the  rate  of  dissipation  of  concentration  fluctuations.  This  surface 
layer  dissipation  is  simulated  in  the  model  by  an  image  sink  for  variance 
which  has  a strength  equal  to  a fraction  a of  the  variance  source  strength. 

It  is  useful  to  consider  the  vertical  source  displacement  and  the  image 
sink  terms  separately.  For  a source  which  is  far  above  the  surface,  where 
neither  the  mean  concentration  image  reflection  source  nor  the  variance  image 
dissipation  sink  have  any  effect,  it  is  easy  to  show  from  Wilson  (1982)  that 


i 

oo 


exp 


exp 


r(z-h)^ 


(3-13) 


where  i is  the  intensity  in  homogeneous  turbulence,  and  i is  the  intensity 
in  the  absence  of  image  sources  and  sinks.  Note  that  i®  includes  the 
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apparent  vertical  displacement  of  the  variance  source 

(3-14) 

with  (p^  = 0.706.  Wilson  (1982)  gives  the  relationship  between  the  intensity 
i at  any  location  with  the  value  i^  with  no  surface  dissipation  as 
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1 - a exp 

2h  z 
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At  the  ground  surface,  z = 0 and  this  reduces  to 
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M 1 0, 5 
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(3-16) 
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= 0.316,  for  a = 0.6 

i ndependen t 

of  source  height.  On  the  centerline  of  a ground  level 

source 
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0 and  (3-13)  reduces  to 

i 

CD 

h=0 

= 1 exp  [- 

(3-17) 

= 0.883  i,  for  (1)  = 0.706 


Combining  the  surface  dissipation  and  wind  shear  source  displacement  effects 
in  (3-16)  and  (3-17)  the  ground  level  intensity  from  a ground  level  source  is 

A,  A 

0.28  i,  where  i is  the  plume  axis  intensity  in  homogeneous  turbulence,  from 
(3-1)  or  (3-10). 

Finally,  it  is  interesting  to  compare  this  simple  model  with  the  results 
of  Sykes,  Lewellen  and  Parker  (1984)  who  solved  the  complete  differential 
equation  for  the  variance  field,  using  closure  models  for  the  various  length 
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2 

and  velocity  scales.  They  found  that  the  boundary  condition  of  5c'  /5z  = 0 
rapidly  developed  into  profiles  with  steep  gradients  close  to  the  surface. 
Their  profiles  were  almost  identical  to  the  present  model  for  both  elevated 
and  ground  level  sources.  (Compare  Figures  7 and  8 of  Sykes,  Lewellen  and 
Parker  (1984)  to  Figures  2 and  4 of  Wilson,  Fackrell  and  Robins  (1982)). 


Intensity  of  Conditionally  Averaged  Fluctuations 

Wilson  (1982)  and  Wilson,  Robins  and  Fackrell  (1985)  develop  a model  for 
the  conditionally  averaged  "plume"  intensity,  in  which  the  intermittent 
periods  of  zero  concentration  are  removed.  The  argument  on  which  this  model 
is  based  is  that  the  total  fluctuation  intensity  can  be  used  to  simulate  the 
conditional  intensity  by  a judicious  choice  of  the  effective  turbulence  scale 
A or  in  (3-1)  or  (3-10).  The  experimental  data  shows  that  for  a ground 

level  source  y = 1.0,  on  the  plume  centerline  from  which  i^  = i there. 

For  the  model  in  (3-1)  and  (3-2),  Wilson,  Robins  and  Fackrell  (1985) 
show  that  the  appropriate  choice  of  turbulence  scale  to  simulate  conditional 
averaging  is  A = 0,  from  which  we  obtain  = 0.12  and  the  plume  intensity 


0.34  X 


p 


(0.12 


(3-18) 


For  the  inertial  subrange  model  in  (3-10),  and  (3-11),  a ground  level  source 


with  h = 0 has  z^^^  = 0.5  XH.  Using  the  scale  relationship  A^^^  = 0.6  z^^^ 

derived  from  the  wind  tunnel  data,  we  find  that  = 0.3  X for  this  case 

Using  this  scale  in  (3-10)  and  (3-11),  the  inertial  subrange  model  predicts 
the  conditional  average  intensity 
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1.37 
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(3-19) 


where  the  virtual  origin  for  the  conditional  intensity  is  at 


32 


1/3  . 2/3 

\ = 0.54  \ (-)  (3-20) 

sp 

There  are  significant  differences  between  the  conditional  "plume"  intensities 
predicted  by  the  empirical  model  in  (3-18)  and  by  the  inertial  subrange  model 
in  (3-19).  The  most  important  difference  is  that  (3-19)  and  (3-20)  predict  a 
weak  dependence  of  conditional  intensity  on  source  size  for  moderate  travel 
distances.  The  model  in  (3-18)  was  purposely  adjusted  to  have  no  source  size 
dependence  for  conditional  intensity  i^. 

Comparison  with  the  experimental  data  in  Figure  (3-3)  supports  the  weak 
dependence  of  i^  on  source  size.  The  inertial  subrange  model  also  shows  a 
much  more  rapid  increase  in  conditional  intensity  from  its  initial  value  of 
zero  at  the  source.  This  is  a considerable  improvement  over  the  Wilson, 
Robins  and  Fackrell  (1985)  prediction  using  (3-18)  and  shown  in  Figure  (3-4) 
for  comparison. 

Effective  Source  Size  for  a Vertical  Jet 

What  value  should  we  use  for  the  source  size  of  a high  velocity  jet  that 
is  bent  over  by  the  wind?  Clearly,  atmospheric  turbulence  will  not  produce 
significant  meandering  or  entrainment  until  the  jet  is  slowed  by  its  own 
self-generated  turbulent  entrainment  to  a speed  where  it  can  be  bent  over  by 
the  wind.  As  a rough  estimate,  we  will  assume  that  the  jet  acts  like  a 
passive  isokinetic  release  once  its  vertical  velocity  W is  equal  to  the  mean 
windspeed  U^.  At  this  point  the  jet  will  be  bent  over  at  an  angle  of  45°. 

The  momentum  jet  analysis  in  Appendix  C shows  that  this  occurs  when  the 
plume  radius  R is  given  by  (C36).  The  effective  source  diameter  d^^^  for  a 
vertical  jet  of  velocity  W and  density  injected  into  a crosswind  of 
velocity  U and  density  is  from  (C36) 
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(3-21) 


The  diameter  is  the  diameter  of  the  "source"  jet  after  any  compressible 
sudden  expansion  needed  to  adjust  the  exit  pressure  to  atmospheric  pressure. 
This  expansion  process  is  discussed  by  Wilson  (1981b)  p.  19,  who  showed  that 
for  a pipe  exit  pressure  greater  than  about  five  times  the  atmospheric 
pressure  that 
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a_ 


(3-22) 


where  d is  the  actual  exit  diameter  (d  = 2R^  on  Figure  4-1)  and  k is  the 
specific  heat  ratio;  k = 1.4  for  air. 

The  diameter  d^^^  calculated  for  (3-21)  is  used  instead  of  the 
isokinetic  diameter  d in  (3-11)  or  (3-20)  to  compute  the  virtual  origin  terms 


\ and  \ for  concentration  fluctuations.  For  a high  pressure  release  from 
a pipeline  rupture  the  difference  between  d and  d^^^  can  be  large.  Using 
Wilson  (1981b)  for  P^/P^  = 10,  d^  - 2d  in  (3-22),  and  with  typically 

100  or  greater,  from  (3-21)  we  find  d^^^  - 200  d,  which  is  certainly  a 
significant  correction. 


Summary 

A new  model  for  along-wind  development  of  concentration  fluctuation 
intensity  has  been  proposed,  based  on  diffusion  in  the  -5/3  inertial  subrange 
of  the  turbulence  spectrum.  The  model  accounts  explicitly  for  source  size 
effects  rather  than  relying  on  a purely  empirical  formulation.  Although  it 
is  somewhat  less  accurate  than  the  purely  empirical  correlation  in  predicting 
peak  intensities,  it  produces  more  realistic  values  close  to  and  far  from  the 


source. 


CHAPTER  4 


DILUTION  CAUSED  RY  SELF-GENERATED  PLUME  TURBULENCE 

The  high  velocity  jet  from  a pipe  rupture  undergoes  three  different  stages 
of  dilution,  as  shown  in  Figure  4-1.  At  the  pipe  exit  the  flow  is  choked  at 
Mach  number  of  unity.  Immediately  outside  the  exit  the  high  pressure  jet 
expands  to  reach  pressure  equilibrium  with  the  atmosphere.  During  this  rapid 
expansion  phase  the  following  analysis  assumes  that  there  is  no  entrainment  of 
ambient  air  by  the  jet.  Following  the  jet  expansion  the  three  phase  of  dilutior 
are: 

1)  The  high  velocity  jet  entrains  air  by  self  generated  jet  turbulence. 
The  entrained  air  slows  the  jet  until  it  can  bent  over  by  the  wind. 

2)  The  bent  over  jet  continues  to  entrain  air  as  it  rises,  and 

coun ter-rotati ng  vortices  are  induced  by  the  process  of  bending  over. 
This  rise  of  the  bent  over  plume  is  perpendicular  rather  than  paralle 
to  its  axis,  and  this  enhances  the  rate  of  ambient  air  entrainment 
into  the  plume. 

3)  The  plume  reaches  its  final  rise  height,  and  further  dilution 
continues  by  the  process  of  turbulent  diffusion,  caused  by  atmospheri 
turbul ence. 

In  the  following  sections,  we  will  make  estimates  of  the  dilution  caused  by 
self-generated  turbulence  in  the  high  velocity  source  jet  and  the  rising  bent 
over  plume,  and  show  how  they  can  be  superposed  on  dilution  from  atmospheric 
turbul ence. 

Dilution  by  Self-Generated  Turbulence  in  the  Exhaust  Jet 

For  an  incompressible  jet  the  pressure  P at  the  pipe  exit  will  be  in 
equilibrium  with  the  atmospheric  pressure  . The  sudden  expansion  from  R^to  F 
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Figure  4-1  Expansion  and  entrainment  into  a bent-over  buoyant  jet  from  a high  pressure 
pipe  release 
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shown  in  Figure  4-1  will  not  occur  in  incompressible  flow,  so  that  and 

can  be  used  in  the  general  compressible  jet  dilution  equation  (C7), 

which  is 

0.5 

4-  = 1 + 2a  (^)  -Z-  (4-1 

• p R 

m ^s  s 

e 

where  z is  the  travel  distance,  is  the  radius,  and  a is  the  jet  entrainment 
constant.  The  most  suprising  characteri Stic  of  (4-1)  is  that  the  jet  velocity 
has  no  effect  on  the  mass  flow  dilution  ratio. 

A pipeline  rupture  sometimes  produces  a horizontal  jet  which  is  bounded  by 
the  ground  surface  on  its  lower  side,  resulting  in  a half-conical  jet.  However 

because  the  change  in  geometry  affects  both  the  entrainment  in  equation  (C2)  an 

the  momentum  conservation  in  (C4)  in  the  same  way,  there  is  no  change  in  the 
dilution  equation  (C7),  or  (4-1)  for  this  situation,  except  to  transform  the 
vertical  travel  distance  z to  horizontal  distance,  x. 

For  variable  density  jets  the  direct  entrainment  measurements  of  Ricou  and 
Spalding  (1961)  yield  an  entrainment  constant  a = 0.08.  Another  careful  set  of 
measurements  by  Crow  and  Champagne  (1971)  integrated  measured  velocity  profiles 
at  varying  downwind  distances  and  obtained  a = 0.073.  Most  investi gators  repor 
values  of  a that  range  from  0.08  to  0.11.  A value  of  a = 0.08  is  a good 
estimate  for  a jet  in  completely  calm  conditions.  However,  even  a small 
crosswind  velocity  U causes  a significant  increase  in  entrainment,  and  a value 
of  a = 0.10  is  recommended  for  jets  discharging  into  a crosswind. 

After  the  vertical  jet  is  bent  over  by  the  wind  it  will  be  carried  along  i 

the  wind  speed  rather  than  its  vertical  velocity  W.  At  the  same  time,  its 
rate  of  entrainment  will  increase  because  the  rise  velocity  W is  perpendicular  j 
to  the  bent-over  jet  axis,  as  shown  in  Figure  4-1.  We  account  for  this  by  an  || 


enhanced  entrainment  constant  p,  which  is  4 to  8 times  larger  than  a. 
flow  dilution  ratio  for  a bent-over  plume  is  given  by  (C16) 
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The  mass 
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(4-2) 


Because  there  is  no  entrainment  during  the  sudden  expansion,  the  mass  flow  at 

the  pipe  exit,  and  at  the  expanded  source  shown  in  Figure  4-1  are  the  same,  so 
2 2 

*^s^^s^s  ~ *^e^e^e*  (^~2)  niay  be  written  in  terms  of  the  atmospheric  pressure 

source  variables  as 
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1 + p (^) 
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(4-3) 


Appendix  C also  develops  equations  for  the  mass  flow  dilution  factors  for 

an  underexpanded  high  pressure  release  from  a pipeline  rupture  or  pressure 

vessel  vent.  For  the  common  situation  where  the  exit  pressure  is  10  times  or 

more  higher  than  atmospheric  pressure  P , the  mass  dilution  ratios  are  given  by 

a 

(C22)  for  a jet  and  (C29)  for  a bent-over  plume. 

Dilution  Factors  for  Contaminant  Concentrations 

The  concentration  of  a toxic  or  flammable  component  in  the  gas  released 
from  a rupture  or  exhaust  vent  is  usually  expressed  as  a volume  fraction  or 
mass  fraction  f^  at  the  pipe  exit.  These  volume  and  mass  fractions  at  the  exit 
are,  for  mixtures  of  ideal  gases,  equal  to  the  values  and  f^  inside  the 
pi  pel ine. 

As  the  gas  in  the  jet  from  the  pipe  exit  entrains  ambient  air,  its  average 
molecular  weight  changes  from  that  of  the  pipe  gas  to  that  of  air  A^.  At  the 
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same  time  the  mass  and  volume  fractions  of  contaminant  decreases  from  f and  C 

e e 

to  f and  C as  exhaust  gas  is  diluted.  The  concentration  dilution  factors  are 
related  to  the  mass  flow  dilution  factor  in  Appendix  C.  If  we  are  far 

enough  from  the  pipe  exit  so  that  entrainment  has  diluted  the  exhaust  gas  to  an 
average  molecular  weight  close  to  air,  the  mass  fraction  dilution  factor  is 
given  by  (C30) 


f 

e = m 


and  the  volume  fraction  dilution  factor  by  (C33) 


(4-4 


(4-5 


The  mass  flow  ratios  may  be  computed  from  (4-1)  or  (4-2)  for  vertical  jets  or 
bent-over  plumes.  In  Appendix  C,  specific  equations  are  derived  for  the  volume 
fraction  dilution  C^/C  from  a high  pressure  compressible  gas  released  as  a 
vertical  jet  (C35),  and  a bent-over  plume  (C34). 


Transition  from  Vertical  Jet  to  Bent-Over  Plume  Behavior 

The  fundamental  problem  with  transition  from  vertical  jet  dilution  in  (4-1 
to  bent-over  plume  dilution  in  (4-2)  is  that  the  so-called  entrainment 
"constants"  a and  p are  not  really  constants.  In  a real  plume  (see  Figure 
4-1)  the  entrainment  "constant"  increases  smoothly  from  a - 0.08  near  the  pipe 
exit  to  p ==  0.6  as  the  plume  bends  over  and  changes  its  entrainment  character- 
istics. Because  (4-1)  and  (4-2)  were  derived  from  the  assumption  of  constant 
values  of  a and  p,  it  is  a mathematical  form  of  self-delusion  to  insert  variab 
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functions  for  a and  p to  force  the  predictions  of  (4-1)  and  (4-2)  to  blend 
smoothly  during  transition. 

At  what  point  should  we  switch  from  the  jet  dilution  of  (4-1)  to  the 
bent-over  plume  dilution  of  (4-2)?  Because  our  interest  is  in  predicting  ground 
level  concentrations  of  toxic  gas,  we  need  only  concern  ourselves  with  the 
bent-over  plume.  It  is  reasonable  to  assume  that  the  increased  entrainment 
constant  p,  and  dilution  predicted  by  (4-2),  will  apply  once  the  plume  has  been 
bent  over  to  a 45°  angle.  As  shown  in  Appendix  C,  this  will  occur  when 
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transi tion 

which  is  (C39)  with  z = Ah,  the  plume  rise.  The  jet  dilution  equation  (4-1) 
will  apply  for  plume  rises  smaller  than  the  transition  value  in  (4-6),  while  the 
bent-over  plume  dilution  of  (4-2)  is  used  for  larger  rises.  It  should  be  kept 
in  mind  that  the  dilution  predicted  by  (4-1)  for  a jet  and  (4-2)  for  a bent-over 
plume  will  not  be  equal  at  the  point  of  transition. 

How  diluted  is  the  plume  once  it  begins  to  bend  over?  This  can  be  estimat- 
ed using  the  transition  condition  (4-6)  in  (4-3),  to  show 


transi ti on 


W 0.5^  2 
h 


(4-7) 


e u — j 

For  high  pressure  ruptures,  > 10  is  a safe  assumption,  and  with  a = 0.1 

and  p = 0.6,  (4-6)  predicts  that  the  mass  flow  dilution  will  always  be  greater 
than  100  as  the  plume  enters  its  bent  over  phase.  It  is  clear  that  self- 
generated turbulence  provides  considerable  dilution  near  the  source. 
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Combined  Dilution  from  Self-Generated  and  Atmospheric  Turbulence 

It  is  essential  that  our  methods  for  estimating  dilution  by  self-generated 
turbulence  produce  a smooth  transition  to  the  dilution  predicted  by  the 
atmospheric  turbulence  dispersion  model.  This  Gaussian  dispersion  equation 
gives  the  plume  centerline  concentration  as 
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(4-8 


in  the  absence  of  ground  surface  reflection. 


The  problem  is  to  relate  this  Gaussian  plume  model  with  spreads  and 


to  the  uniform  concentra tion  profile  of  radius  R used  in  plume  growth  from  self 
generated  turbulence.  Wilson  (1979)  suggests  that  for  the  conical  plume  from 
self-generated  turbulence  that  the  plume  radius  is 


R = 2 a 


sel  f 


(4-‘ 


It  is  easy  to  show  that  this  radius  implies  a uniform  concentration  of  0.5 
across  the  top-hat  plume  profile  to  have  the  same  contaminant  flow  rate  Q as  tl 
Gaussian  in  (4-8).  Using  (CIO)  for  R = R + pAh,  and  (C15)  to  estimate  the 
effective  initial  radius  R^  of  the  bent  over  plume  we  find  from  (4-9) 
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Typically  the  pipe  radius  R^  is  small,  and  - 100  for  high  pressure  sonic 

flow  at  the  exit.  In  this  case  the  initial  radius  R - 5 R which  is  negligib 

0 e 

compared  to  the  pAh/2  term.  The  self-generated  spread  becomes 


pAh 

^self  ^ 
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as  suggested  by  Wilson  (1979).  Then,  assuming  that  the  self-generated 
turbulence  is  completely  uncorrelated  with  atmospheric  turbulence,  the  spreading 
variances  simply  add,  and 


0.5 


sel  f ' 


(4-12) 


where  is  the  combined  atmospheric  and  self-generated  spread. 

How  important  is  the  spread  from  self-generated  turbulence?  For  ground 
level  concentrations  from  high  stacks,  always  negligible  compared  to 

a^.  However,  for  a ground  level  source  the  situation  might  be  different.  For 
the  special  case  where  ^ constant,  the  maximum  ground  level 

concentration  in  a Gaussian  plume  from  a ground  level  source  occurs  when 


(4-13) 


Using  (4-11)  in  (4-12)  and  then  equating  this  to  (4-13)  leads  to 
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at  the  point  of  maximum  ground  level  concentration.  For  p = 0.6  we  find  that 

^sel  f ” ^zc  ^ ^z*  ^ increase,  and  given 

the  high  level  of  uncertainty  in  estimating  in  the  atmosphere,  ^ can  also 

be  neglected  for  ground  level  sources  with  plume  rise. 

Because  the  magnitude  of  this  ground  level  maximum  concentration  is 

proportional  to  the  ratio  ^^c'^'^yc’  self-generated  spread  will  have  a 

negligible  effect  because  it  increases  both  o and  a . From  this,  we  conclude 

y c zc 

that  self-generated  turbulence  within  the  plume  is  unlikely  to  make  any 
significant  contribution  to  the  location  of  maximum  ground  level  concentration. 


or  to  its  magnitude.  In  contrast,  self-generated  spread  is  the  major  factor 
contributing  to  plume  dilution  during  its  initial  vertical  jet  and  bent  over 
phases,  as  shown  by  (4-7). 


CHAPTER  5 


MOMENTUM  AND  BUOYANCY  CONTRIBUTIONS  TO  FINAL  PLUME  RISE 


Accidental  releases  of  toxic  and  flammable  gases  often  occur  as  high  velo- 
city jets  from  pressure  vessels  and  pipelines.  They  differ  from  stack  gas  emis- 
sions in  two  important  ways:  jet  momentum  makes  a significant  contribution  to 

the  trajectory  and  to  the  final  rise  height  of  these  plumes,  and  secondly,  these 
releases  often  occur  close  to  ground  level  and  pass  through  the  varying  wind- 
speed  and  turbulence  of  the  surface  boundary  layer.  These  differences  require 
changes  in  the  way  plume  rise  is  evaluated,  in  particular,  the  way  in  which  the 
effects  of  momentum  and  buoyancy  are  superposed  to  determine  combined  rise. 

Briggs  (1975)  presents  a solution  for  plume  rise  for  the  case  where  the 
entrainment  "constant"  p is  a true  constant  with  downwind  distance.  For  a plume 
with  an  initially  small  diameter  at  the  point  of  release,  the  plume  rise  during 
this  phase,  where  self-generated  turbulence  determines  the  rate  of  entrainment 
of  ambient  air. 


That  is,  that  the  contributions  of  momentum  and  buoyancy  add  as  a sum  of  cubes 
in  the  transitional  rise  phase. 

To  predict  final  plume  rise  Wilson  (1981)  suggested  that  the  sum  of  cubes 
law  be  used  as  an  approximate  way  of  combining  the  final  rises  predicted  for  the 
separate  momentum  and  buoyancy  components.  This  approach  is  only  an  approxima- 
tion because  in  the  actual  combined  momentum- buoyancy  plume  there  is  a single 
downwind  location  that  corresponds  to  final  rise,  while  in  the  approximation 
downwind  distances  to  final  momentum  rise  and  final  buoyancy  rise  are  generally 


1/3 


(5-1) 


-45- 


46 


different.  The  sum  of  cubes  law  in  (5-1)  is  only  strictly  valid  for  rise  com- 
ponents at  the  same  downwind  location,  and  not  for  separate  components  treated 
as  individual  plumes  with  different  downwind  distances  to  final  rise. 

In  this  chapter  we  will  examine  the  errors  associated  with  using  this  sum 
of  cubes  law  as  an  approximation  for  combining  the  separate  momentum  and  buoy- 
ancy contributions  to  final  rise.  Because  the  errors  in  the  approximation 
depend  on  the  theoretical  model  chosen  to  predict  the  point  of  final  rise,  it  is 
worthwhile  to  first  review  some  of  the  available  physical  and  theoretical  models 
for  this  process. 

Factors  Determining  Final  Plume  Rise 

In  a stable  atmosphere,  the  approach  to  final  rise  is  a well  understood 
process.  A stable  atmospheric  temperature  variation  with  height  produces  a 
downward  force  on  the  rising  plume  which  eventually  causes  it  to  level  off  at 
its  final  rise  height.  In  an  unstable  or  neutrally  stable  atmosphere  this 
restoring  force  is  absent,  and  because  momentum  and  buoyancy  fluxes  are  conserv- 
ed, the  plume  should  continue  to  rise  indefinitely  as  it  is  carried  downwind. 
However,  when  the  atmosphere  is  turbulent,  plumes  tend  to  level  out  at  a final 
rise  height  regardless  of  atmospheric  stability.  While  it  is  clear  that  some 
interaction  occurs  between  the  plume  and  the  atmospheric  turbulence  which 
produces  this  final  rise,  there  is  no  general  agreement  as  to  the  physical 
mechanism  by  which  this  occurs.  There  appear  to  be  three  distinctly  different 
schools  of  thought  on  the  process  which  leads  to  final  rise.  These  three 
physical  models  are: 

1.  Atmospheric  turbulence  gradually  modifies  the  entrainment  con- 
stant, causing  p to  increase  with  downwind  distance.  The  plume 
rise  never  actually  terminates,  but  the  very  large  entrainment 
constant  p causes  the  rate  of  rise  to  be  negligibly  small  far 
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downwind.  This  approach  was  used  by  Bosanquet,  Carey  and  Hal  ton 
(1950),  Bosanquet  (1957),  Briggs  (1969)  and  Slawson  and  Csanady 
(1967),  (1971). 

2.  Atmospheric  turbulence  causes  an  "extrai nmen t"  of  material  from 
the  "dynamic"  plume.  This  extrained  fluid  ceases  to  assist  the 
plume  in  rising  as  a coherent  entity.  Extrainment  is  a gradual 
process  which  eventually  completely  depletes  the  plume  of 
"active"  or  "dynamic"  buoyancy  and  momentum  flux,  and  leads  to  a 
leveling  off  at  final  rise.  The  use  of  extrainment  to  deplete 
the  plume  of  buoyancy  and  vertical  momentum  has  been  used  in 
models  proposed  by  Priestly  (1956),  Turner  (1963),  Benjamin 
(1975)  and  Netterville  (1983). 

3.  Atmospheric  turbulence  has  little  effect  on  the  entrainment  pro- 
cess until  some  critical  point  is  reached  where  the  velocity 
scale  of  internal,  sel f-generated  plume  turbulence  becomes  less 
than  the  equivalent  atmospheric  turbulence  velocity  scale.  At 
that  point,  atmospheric  turbulence  "breaks  up"  the  plume  into 
uncorrelated  parcels  which  cease  to  rise  in  an  organized  way. 

This  simple  model  in  which  the  transitional  rise  is  unaffected  by 
atmospheric  turbulence,  and  plume  break-up  occurs  so  rapidly  that 
there  is  a negligible  transition  distance  between  transitional 
and  final  rise,  has  been  used  by  Scorer  (1959)  and  by  Briggs 
(1975). 


Briggs  (1969),  (1975)  presents  a considerable  amount  of  experimental  evi- 
dence to  show  that  buoyant  plumes  follow  the  transitional  rise  law  given  in 
(5-1)  over  a large  portion  of  their  trajectory.  For  this  reason  the  simple 
concept  of  unimpaired  transitional  rise  followed  by  rapid  plume  breakup  is  a 
useful  one  for  predicting  both  trajectory  and  final  rise  height.  In  the 
following  section  a new  model  for  final  rise  based  on  local  velocity  scales 
within  the  inertial  subrange  of  the  atmospheric  turbulent  spectrum  will  be 
developed.  We  will  show  that  this  model  is  equivalent  to  Briggs  approach  of 
defining  a characteri Stic  velocity  scale  by  equating  the  rate  of  dissipation  of 
internal  sel f-generated  plume  turbulence  to  that  of  the  atmospheric  turbulence. 


Plume  Breakup  Criteria  From  Local  Spectral  Scaling 


Before  we  propose  a new  criteria  for  the  point  of  plume  breakup  and  final 
rise  let's  look  at  the  existing  models  of  Scorer  and  Briggs.  Scorer  (1959)  made 
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the  simple  assumption  that  plume  breakup  occurs  at  the  point  where  the  RMS 
atmospheric  turbulence  intensity  u‘  was  equal  to  the  average  rise  velocity  W o 
the  plume.  In  the  surface  layer  of  the  atmosphere  where  the  RMS  intensity  u' 
remains  constant  with  height.  Scorer's  model  leads  to  a final  buoyant  rise  giv 
by 


Ah 


Bf 


Scorer's  model  is  unrealistic  because  it  assumes  that  all  of  the  atmospheric 
turbulence  energy  contributes  to  plume  breakup.  For  a real  plume,  the  large 
scale  atmospheric  turbulence  acts  only  to  cause  plume  meandering  but  does  not 
break  up  the  organized  motions  that  contribute  to  plume  rise.  The  atmospheric 
turbulence  eddies  which  break  up  the  plume  must  have  a scale  about  the  same  si 
as  the  plume  itself,  characterized  by  its  radius  R. 

Briggs  (1975)  accounted  for  the  scales  of  atmospheric  and  plume  turbulen( 
by  comparing  the  plume's  inertial  subrange  turbulence  with  that  of  the  atmos- 
phere. By  considering  only  the  inertial  subrange,  large  scale  motions  that 
contribute  to  meandering  were  effectively  excluded.  The  turbulence  spectrum 
the  inertial  subrange  is  characteri zed  by  the  rate  of  energy  transfer  which  i! 
equal  to  e,  the  turbulence  dissipation.  Using  this,  Briggs  set  the  point  of  I 
plume  breakup  at  the  location  where  the  dissipation  rate  of  internal  self-  I 
generated  turbulence  within  the  plume,  F equal  to  the  atmospheric  dissipation  I 
e.  Briggs'  criteria  for  equality  of  dissipation  rates  at  final  rise  is  I 

discussed  in  detail  in  Appendix  D where  the  final  rise  of  a purely  buoyant  pi  I 

3 ' 

is  given  in  (D17)  and  (028)  which  for  e = u*/kz  and  the  friction  velocity  u* 
proportional  to  the  windspeed  U,  leads  to 
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this  is  identical  to  (5-2)  for  a ground  level  release  into  the  roughness 
generated  turbulence  of  the  surface  layer.  At  the  other  extreme, 
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for  high  stacks  with  small  plume  rise,  in  which  case  the  atmospheric  dissipation 

3 

has  an  approximately  constant  value  of  e - u*/i<h^  through  the  plume  layer.  It 
is  interesting  to  note  that  neither  of  these  final  rise  predictions  correspond 
to 

i:3/5 

Ahpf  . ^ (5-5) 

developed  by  Briggs  (1969)  and  currently  recommended  by  Alberta  Environment  for 
stack  design. 

The  major  difficulty  with  Briggs  approach  is  a conceptual  one.  It  is  dif- 
ficult to  visualize  the  physical  processes  involved  in  breakup  when  the  turbu- 
lence energy  dissipation  rate  is  used  as  the  governing  parameter.  The  basic 
concept  of  rapid  plume  breakup  begins  by  assuming  that  the  plume  radius  is 
initially  small,  and,  because  there  is  little  turbulence  energy  in  atmospheric 
scales  of  this  size  the  plume  growth  is  dominated  by  self-generated  turbulence. 
As  the  plume  rises  and  entrains  ambient  air  its  rise  velocity  and  its  internal 
turbulence  velocity  scale  both  decrease,  and  its  length  scale  of  turbulence 
increases  in  proportion  to  R,  the  plume  radius.  After  the  plume  has  grown  to 
sufficiently  large  size,  atmospheric  eddies  of  scale  R have  the  same  turbulence 
energy  as  self-generated  plume  turbulence.  However,  because  these  entrained 
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eddies  have  motions  that  are  uncorrelated  with  the  plume  turbulence  they  cause 
the  plume  eddies  to  lose  their  correlation,  and  the  organized  rising  motion  of 
the  plume  breaks  up. 

A limiting  condition  for  plume  breakup  that  is  physically  more  satisfying 


than  Briggs'  equality  of  dissipation  rates  is  to  equate  the  turbulence  energy 


2 2 

w'  of  self  generated  plume  turbulence  to  the  turbulence  energy  w'.  of  that  par 
P ^ 


of  the  atmospheric  spectrum  that  has  the  same  length  scale  of  turbulence  as 


the  plume.  The  subscript  "k"  denotes  the  energy  and  scale  at  a particular 
wavenumber  k in  the  spectrum.  (The  wavenumber  k is  a normalized  fluctuation 
frequency  which  is  inversely  proportional  to  the  wavelength  of  the  turbulence’ 
Our  assumption  is  that  plume  breakup  occurs  when 
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Using  the  principle  of  local  scaling  of  the  turbulence  spectrum  developed  in 
Appendix  A as  (A13)  and  (A14)  we  can  write  the  RMS  velocities  at  turbulence 


scale  as 
k 


"k  = 


L 
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where  A^  is  a constant  and  z is  the  turbulence  dissipation  rate.  In  the  same 

way,  for  the  internally  generated  plume  turbulence  with  its  dissipation  rate? 
and  scale  Z 

P 


(5-8 


The  relevant  part  of  the  atmospheric  turbulence  spectrum  that  influences  break 
will  have  the  same  scale  as  the  plume  turbulence.  Setting  Z,  = Z - R,  the 

K p 
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plume  radius,  in  (5-7)  and  (5-8)  and  inserting  these  in  (5-6)  results  in  a plume 
breakup  condition  of 
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(5-9) 


This  is  exactly  the  same  as  Briggs'  dissipation  condition,  except  that  it  is 

less  sensitive  to  a mismatch  in  dissipation  because  for  equal  energy, 

dissipation  e appears  to  the  2/3  power. 

The  use  of  local  spectral  scaling  for  plume  dissipation  e is  discussed  in 

— 3 

Appendix  D where  it  is  shown  that  Briggs'  empirical  constant  g in  e = g W /Ah 
can  be  derived  from  other  known  constants  related  to  the  turbulence  spectrum. 
The  good  agreement  between  the  derived  value  of  n = 1.33  and  Briggs'  empirical 
value  of  g = 1.5  is  further  evidence  for  the  validity  of  local  spectral  scaling 
to  define  plume  turbulence  velocity  and  length  scales. 


Sum  of  Cubes  Approximation  to  Final  Rise 

The  final  rise  height  of  a plume  with  combined  momentum  and  buoyancy  can  be 
determined  exactly  by  using  the  sum  of  cubes  law  in  (5-1).  The  contributions 
from  momentum  and  buoyancy  are  evaluated  from  (D31)  and  (D32)  in  Appendix  D, 
with  the  time  of  travel  set  at  t^,  the  final  rise  time  from  (D25)  or  (D28).  The 
theory  for  determining  the  final  rise  time  using  Briggs  dissipation  criteria  for 
plume  breakup  is  developed  in  Appendix  D. 

With  an  exact  expression  for  final  rise  in  (5-1),  why  bother  developing  an 
approximation?  The  answer  lies  in  the  functional  form  for  the  final  rise  time 
t^  given  in  its  general  form  in  (Dll)  and  for  specific  vertical  variations  of 
atmospheric  dissipation  in  (D25)  and  (D28).  These  functional  forms  for  the 
final  rise  time  must  be  solved  by  iteration,  which  precludes  any  simple  closed 
form  solution  for  the  relative  contributions  of  buoyancy  and  the  momentum. 
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It  is  simpler  and  physically  more  sati  si  tying  to  compute  the  final  buoyanc 
and  momentum  rises  separately  and  combine  them  in  some  way  to  determine  the 
final  rise.  Because  the  separate  final  rises  occur  at  different  travel  times, 
sum  of  cubes  of  separate  final  rises  similar  to  (5-1)  is  only  an  approximation 
to  the  actual  final  rise.  The  nature  of  this  approximation  is  easiest  to  see  t 
looking  at  the  travel  time  to  final  rise  of  a purely  buoyant  plume,  a momentum 
jet,  and  a combined  momentum-buoyant  plume.  In  the  combined  plume  both  momenti 
and  buoyancy  contribute  to  the  production  of  self-generated  turbulence.  The 
turbulence  energy  and  dissipation  rate  within  this  combined  plume  will  always  I 
larger  than  for  a purely  buoyant  plume  or  a pure  momentum  jet.  The  larger 
turbulence  energy  in  the  combined  plume  allows  it  to  resist  breakup  by  atmos- 
pheric turbulence  for  a greater  time  of  travel.  From  this  we  see  that  the  tim 
of  travel  to  final  rise  of  a combined  plume  will  always  be  larger  than  the 
momentum  jet  component  or  the  buoyant  plume  component. 

However,  from  a practical  point  of  view,  the  relevant  question  is  how 
accurate  a sum  of  cubes  of  separate  final  rises  is  as  an  approximation  to  the 
combined  final  rise  of  (5-1).  This  approximation  is  discussed  in  detail  in 
Appendix  D where  the  approximate  sum  of  cubes  of  separate  final  rises  is  wri tt 
as 


(5-1 


where  the  overbar  " indicates  the  separate  final  rises  of  the  momentum  jet 
and  buoyant  plume  components,  and  Ah^  is  the  sum  of  cubes  of  separate  rises. 

The  separate  momentum  and  buoyancy  final  rises  Ah^^^  and  Ah^^^  are  derived 
Appendix  D as  (016)  and  (017)  with  (021)  for  constant  dissipation  e through  tl^ 
plume  layer.  For  a ground  level  release  into  roughness  generated  turbulence  ' 
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Figure  5-1  Error  induced  by  using  sum  of  cubes  law  for  separate  momentum  and  buoyant 
final  rises  instead  of  combined  final  rise.  Theory  from  (D25),  (D35)  and 
(D28),  (D36). 


momentum  and  buoyancy  rises  are  given  by  (D16)  and  (D17)  with  (022)  for  the 
dissipation  e. 

The  ratio  Ah^/Ah^  of  the  exact  and  approximate  final  rises  in  (5-1)  at  t^ 
and  in  (5-12)  for  separate  final  rise  times  tj^^  and  t^^  is  given  in  Appendix  0 
by  (035)  and  (036).  This  ratio,  which  represents  the  error  in  the  "sum  of  cube 
of  separate  rises"  approximation,  has  the  useful  character! Stic  that  it  can  be 
made  a function  only  of  the  ratio  of  buoyant  to  momentum  rises  Ah^^/Ahj^^.  In 
this  form  it  is  independent  of  the  entrainment  constant  p and  the  sel f-generate 
turbulence  dissipation  rate  constant  g,  both  of  which  must  be  estimated  from 
experiment. 

The  ratio  of  the  exact  to  approximate  sum  of  cubes  plume  rises  from  (D25) 
(D35)  and  (D28),  (D36)  are  shown  in  Figure  5-1.  The  largest  errors  occur  when 
the  final  momentum  and  buoyancy  rises  are  approximately  equal.  For  high  stack 
where  the  atmospheric  dissipation  is  constant  through  the  plume  layer,  the  err 
is  negligible  at  all  buoyancy  to  momentum  rise  ratios.  For  the  extreme  case  o 
a ground  level  release  into  roughness  generated  surface  layer  turbulence,  the 
atmospheric  dissipation  is  inversely  proportional  to  the  height  above  ground, 
see  (D19),  and  errors  of  about  10%  result  from  the  sum  of  cubes  approximation. 
However,  the  most  important  result  shown  in  Figure  5-1  is  that  the  approximati 
using  the  sum  of  cubes  of  separate  final  rises  will  usually  underestimate  the 
actual  plume  rise.  This  will  cause  diffusion  models  to  slightly  overestimate 
ground  level  concentrations,  so  that  the  separate  rises  approximation  will  err 
slightly  on  the  conservative  side.  H 
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TOTAL  AND  INSTANTANEOUS  PLUME  SPREAD  DETERMINED  DIRECTLY 
FROM  TURBULENCE  PARAMETERS 

To  estimate  the  concentration  fluctuation  intensity  for  varying 
sampling  times,  and  at  different  downwind  distances  from  the  source,  it  is 
necessary  to  specify  the  instantaneous  plume  widths  a-jy  and  a-j^  and  the 
time  averaged  "total"  widths  ay  and  a^.  To  be  useful  for  a wide  range 
of  atmospheric  conditions,  these  plume  widths  should  be  in  terms  of  turb- 
ulence length  and  velocity  scales  that  can  be  estimated  or  directly 
measured.  Simple  functional  forms  for  a and  a-j  that  are  consistent  with 
the  required  asymptotic  behavior  near  to  and  far  from  the  source  will  be 
developed  in  the  following  sections. 


Total  Plume  Spread 

To  simplify  the  presentation  only  the  crosswind  width  oy  will  be 
considered  in  detail,  and  the  equivalent  results  for  the  vertical  spread 
a^  written  by  analogy.  The  total  time  averaged  plume  width  used  to 
estimate  the  mean  concentration  c is  made  up  of  an  instantaneous  component 
a^.^,  and  a meandering  component  a^^,  which  are  assumed  to  be  uncorrelated 
and  so  have  additive  variances 
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A great  deal  of  field  and  laboratory  measurements  and  analysis  exist  to 
estimate  this  total  spread.  The  measurements  suggest  that  over  a very  wide 
range  of  downwind  distances  the  plume  spread  follows  a power  law  of  the 
form 


O 2 = n 2 , 2 

y y s 


(A2) 


-55- 
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The  source  size  is  usually  regarded  as  negligible.  Because  of  the 
ease  with  which  power  law  expressions  can  be  used  in  diffusion  analysis,  we 
will  propose  a model  based  loosely  on  Taylor's  (1921)  statistical  theory 
that  preserves  this  power  law  form.  Assuming  that  the  source  size  as 
simply  adds  in  the  variance  we  find,  see  Csanady  (1973)  p.  55-56 

a ^ = (-^)  + a ^ for  x <<  L (A3) 

y s V ^ ^ 

and  ^2  99 

oy  = 2(-^)  (L^x  - Ly^)  + for  x » L (A4) 

Uc 

where  v'  is  the  RMS  turbulence  velocity,  1)^  is  the  mean  advection  speed 
and  Ly  is  the  Lagrangian  integral  length  scale  of  velocity  component  v. 

A simple  interpolation  equation  between  the  near-source  and  far  downwind 
limits  is  given  by  Snyder  (1984)  as 

The  advection  velocity  U^-  is  related  to  the  travel  time  t by  = x/t. 

In  shear  flow  the  advection  velocity  will  in  general  be  different  than  the 
local  velocity  at  source  height.  Some  averaging  process  must  be  used  to 
define  a value  of  that  satisfies  mass  conservation  at  each  plume  cross 
secti on. 

The  Lagrangian  length  scale  Ly  is  usually  estimated  from  the 
Eulerian  integral  scale  of  turbulence  Ay.  This  Eulerian  scale  is  the 
along-wind  correlation  of  the  crosswind  turbulence  velocity  derived  from 
the  time  series  of  velocity  from  a single  wind  sensor,  and  the  assumption 
of  frozen  turbulence.  Corsin  (1963)  suggests  a relation  of  the  form 
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(A6) 


Pasquill  (1974)  pp.  89-92  reviews  the  available  measurements  in  the  atmos- 
pheric boundary  layer  and  finds  values  of  ranging  from  0,35  to  0.80, 
from  which  he  recommends  Ax  = 0.68  0.2.  In  shear-free  grid  turbulence 

Snyder  and  Lumley  (1971)  found  Ax  = 1.0  by  direct  measurement. 


Power  Law  Approximations 

To  interpolate  between  (A3)  and  (A4)  and  at  the  same  time  maintain  a 
power  law  functional  form  like  (A2),  we  assume 


where  from  (A3)  and  (A7),  the  exponent  ay  must  lie  somewhere  between  1.0 
near  the  source  and  0.5  far  downwind.  By  comparing  (A2)  and  (A7)  the 
dimensional  constant  Dy  is 

Dy  = (A8) 

Uc 

Vertical  plume  spread  oy  is  found  in  the  same  way,  replacing  Ay  and  v' 
wi  th  Ay^  and  w' . 

One  remarkable  feature  of  these  power  law  equations  is  the  wide  range 
of  downwind  distances  over  which  they  apply.  Pasquill  (1974)  p.  228 
reviews  some  experiments  that  show  that  crosswind  spread  follows  a power 
law  with  ay  = 0.8  for  distances  up  to  a few  hundred  kilometers.  This 
failure  to  approach  a^  = 0.5  for  distances  larger  than  about  3L^  may  be 
caused  by  crosswind  shear  as  wind  direction  changes  with  height. 

Because  turbulence  scales  are  much  smaller  in  the  vertical,  does 

have  the  expected  a^  - 0.5  for  x > 3L^,  see  Pasquill  (1974)  p.  202.  The 
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process  of  transition  to  this  asymptotic  limit  is  very  gradual,  as  the 

0 5 

growth  varies  from  a « x near  the  source  to  « x far  downwind.  A 
power  law  representation  with  a^  - 0.7  is  often  a reasonable  approximation 
for  distances  from  100  m to  5000  m. 

Instantaneous  Plume  Width 

The  instantaneous  plume  spread  a^y  is  needed  to  estimate  the  effect 
of  meandering  on  both  mean  and  fluctuating  concentrati ons.  Also,  because 
the  spread  of  an  ensemble  of  instantaneous  puffs  about  their  mass  centers 
will  be  the  same  as  the  instantaneous  rather  than  total  spread,  it  is  a-jy 
that  is  used  to  characterize  puff  dispersion. 

The  essential  difference  between  instantaneous  and  total  plume  growth 
is  that  different  ranges  of  eddy  sizes  influence  total  and  instantaneous 
spread.  Neither  a-jy  or  cry  are  significantly  changed  by  eddies  much 
smaller  than  the  instantaneous  spread  a^y.  Motions  of  these  small  eddies 
are  uncorrelated  over  distances  comparable  to  the  instantaneous  width. 

These  uncorrelated  motions  make  little  contribution  to  the  growth  of  a-jy 
or  ay  and  act  mainly  to  smear  out  small  scale  inhomogenieties  within  the 
instantaneous  plume.  Pasquill  (1974)  p.  127  demonstrates  this  effect  by 
spectral  analysis  of  contributions  to  the  total  spread  ay. 

Large  eddies  that  cause  the  instantaneous  plume  to  meander  back  and 
forth  over  a fixed  receptor  make  a contribution  to  total  dispersion  ay 
when  the  sampling  time  is  sufficiently  long  to  include  these  large  time 
scale  motions.  In  contrast,  the  instantaneous  plume  spread  a-jy  repre- 
sents the  process  of  relative  diffusion  about  the  instantaneous  mass  center 
of  a plume  cross  section.  This  spread  is  determined  by  ensemble  averaging 
a large  number  of  instantaneous  plume  realizations  which  have  each  been 
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shifted  to  have  a common  center  of  mass.  Because  the  effect  of  meandering 
by  large  eddies  is  removed  by  considering  dispersion  relative  to  the 
instantaneous  mass  center,  the  growth  of  o^y  is  influenced  only  by  eddies 
about  the  same  size  as  cr-jy.  Large  eddies  act  only  to  push  the  entire 
plume  around,  while  the  uncorrelated  motions  of  small  eddies  act  mostly  to 
reduce  inhomogen iety  within  the  instantaneous  plume. 

In  terms  of  the  spectrum  of  turbulent  velocity  fluctuations  E2(l<)  of 
v'2  the  process  of  "total"  dispersion  that  produces  Oy  is  like  a 
lowpass  filter  that  includes  all  eddies  larger  than  o]y  in  the  diffusion 
process.  On  the  other  hand,  the  relative  dispersion  that  controls  is 
equivalent  to  a bandpass  filter  that  includes  only  narrow  ranges  of  eddy 
sizes  near  the  plume  size  a-jy.  In  the  following  analysis  we  will  exploit 
the  fact  that  only  the  velocity  fluctuations  near  a-jy  need  to  be 
considered  when  estimating  the  instantaneous  spread. 

It  has  been  observed  that  the  instantaneous  spread  a-jy  is  much  more 
sensitive  than  ay  to  the  initial  release  conditions.  Because  of  this,  it 
is  essential  to  consider  the  effect  of  initial  source  size  in  developing 
models  for  instantaneous  width.  We  will  do  this  by  expressing  the  initial 
source  size  as  an  equivalent  virtual  spread 

Estimates  of  Instantaneous  Spread 

A quantitative  estimate  of  instantaneous  plume  spread  in  relative 
diffusion  was  made  by  Batchelor  (1952)  and  later  expanded  by  Batchelor  and 
Townsend  (1956)  using  dimensional  reasoning.  Their  results  apply  to 
turbulence  which  has  an  inertial  subrange  where  E2(k)  falls  off  like 
k"^/^  and  can  be  charac teri zed  by  a single  variable,  the  kinetic  energy 
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dissipation  rate  e.  Following  Pasquill's  (1974)  pp.  141-142  review  of 
their  results,  and  using  large  eddy  scaling  of  dissipation  for  which 

3 

e V /A^,  they  found, 

+ os^  for  x«Ay  (A9) 

Uc  Ay 

very  close  to  the  source.  The  direct  dependence  of  the  growth  rate  on 

Os/Ay  is  due  to  the  source  size  setting  the  initial  width  of  the 

2 

spectral  bandpass  filter  that  determines  the  fraction  of  the  v'  spectrum 
of  that  contributes  to  the  instantaneous  plume  spread.  At  intermediate 
distances  Batchelor  found,  for  negligible  source  size 


a.  2 = AJxL)V 


(AlO) 


This  is  the  accelerated  growth  phase  with  o]  « x ^/2  was  observed 

by  Richardson  (1926).  The  value  of  the  constant  A3  can  be  set  using  the 
results  of  Smith  (1968)  who  applied  a Lagrangian  analysis  of  conditioned 
particle  diffusion  to  show  that  at  intermediate  distances 


a.  ‘ 

ly 


2 3 

2 [iL] 

3 'U  L 

C V 


X « L 


(All) 


Using  (A6)  to  establish  a relation  between  Eulerian  and  Lagrangian  scales 
allows  (All)  to  be  written  in  the  same  form  as  (AlO)  with 


(A12) 


The  accelerated  growth  phase  must  terminate  at  some  distance  because  a-jy 
is  always  less  than  Oy,  which  is  proportional  to  x^y  with  ay  = 1.0 
near  the  source  and  0.5  far  downwind. 
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Spectral  Bandpass  Model  for  Relative  Diffusion 

In  this  section  we  will  propose  a model  for  predicting  the  effect  of 
source  size  on  instantaneous  spread  and  use  it  to  estimate  the  growth  of 
the  instantaneous  plume  in  terms  of  the  turbulence  scales  and  velocities. 

An  added  attraction  of  this  model  will  be  the  physical  insight  it  provides 
into  the  fundamental  differences  between  total  and  i nstan taneous  dispersion 
embodied  in  the  spreads  ay  and  a-jy. 

The  essential  character! Stic  of  relative  diffusion  about  an  instant- 
aneous mass  center  is  that  it  is  the  equivalent  to  imposing  a bandpass 
filter  on  the  turbulent  velocity  spectrum.  Because  of  this  it  is  only 
necessary  to  describe  a relatively  narrow  region  of  the  spectrum  that  is 
active  in  causing  instantaneous  spread.  This  frequency  band  will  be 

^ ~ \ where  the  wavenumber  is  the  inverse  wavelength  at  frequency 

n,  and  is  defined  as 

k = = 2n 

wavelength 

As  grows  larger  the  frequency  band  broadens  and  moves  to  lower  fre- 
quencies until  eventually  the  instantaneous  plume  samples  the  same 
frequencies  as  the  total  spread,  and  is  the  same  as  a^.  In  our  analy- 
sis we  assume  that  the  instantaneous  plume  is  small  enough  that  it  is 
influenced  only  by  eddies  in  the  inertial  subrange.  There  is  no  signifi- 
cant turbulent  dissipation  or  production  of  v'^  in  this  part  of  the  spec- 
trum, so  the  energy  at  each  frequency  depends  only  by  the  rate  at  which 
energy  is  transferred  from  larger  to  smaller  scales.  This  transfer  rate 
must  equal  the  rate  of  dissipation  e which  ultimately  consumes  all  turbu- 
lence kinetic  energy  after  it  passes  through  the  inertial  subrange.  This 


dissipation  can  be  estimated  from  the  turbulence  kinetic  energy  v'^  and 
the  large  eddy  time  scale  t'  = Ay/v' 
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or 


.3 


e = A 


V 


2 * 

where  v'  = (v  ) is  the  RMS  turbulence  velocity,  the  integral  length 


scale  and  A4  is  a constant  of  order  unity. 

If  the  only  controlling  factor  on  the  spectrum  is  the  rate  of  energy 


where  is  the  length  scale  of  eddies  at  wavenumber  k.  Equating  (A13) 
to  (A14)  gives  an  estimate  of  the  RMS  turbulence  at  scale  3S 


This  scaling  has  been  suggested  by  Sykes,  Lewellen  and  Parker  (1984),  and 
will  determine  the  amount  of  RMS  velocity  that  is  effective  in  causing 
instantaneous  plume  spread.  Before  doing  this,  it  is  interesting  to  see 
what  spectrum  shape  is  implied  by  (A15).  Assuming  the  local  eddy  scale 
at  wavenumber  k to  be 


2 

transfer,  then  the  energy  Vj^'  at  any  wavenumber  k must  also  scale  in  the 


same  way  as  (A13),  so  that 


.3 


(A14) 


(A15) 


V 


(A16) 
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where  A5  is  a constant  of  order  unity.  Then  writing  the  definition  of 
spectral  density  E2(k)  of  v'^ 

A(v  '2)  = E-(k)  Ak 


or 


E (k)  = - 
2 dk 


(A17) 


Using  (A16)  in  (A15)  and  evaluating  the  spectrum  function  from  (A17)  leads 
to 


. -2/3  2 -2/3  -5/3 

E2(k)=|A5  V'  k 


(A18) 


Using  (A13)  to  write  v'^  in  terms  of  dissipation  e 


E^{k)  = I (A4A5) 


-2/3  2/3  -5/3 


(A19) 


which  shows  that  the  local  velocity  scaling  of  the  spectrum  assumed  in 
(A14)  leads  to  the  expected  form  for  the  inertial  subrange.  The  value  of 
the  lead  constant  in  (A19)  is  about  0.74,  see  Hinze  (1975)  p.  255.  Using 
the  usual  value  of  A4  = 0.8  we  find  that  A5  = 1.07,  close  to  unity. 


Effective  Eddy  Oiffusivity  for  gjy 

To  model  the  growth  rate  of  the  instantaneous  plume  we  propose  an 

effective  eddy  diffusivity  K^. . The  instantaneous  spread  is  related  to  the 

eddy  diffusivity  by  ^ 

da. 

— = 2 Kiy  (A20) 

dt 

see  Csanady  (1973)  p.  63.  Because  the  instantaneous  plume  responds  only  to 
turbulence  scales  in  the  narrow  spectral  bandpass  near  a-jy,  we  assume  a 
turbulent  diffusivity 
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K.  = A^v.  'a. 
iy  6 1 ly 


for  o.  < A 

ly  V 


(A21) 


where  v-j ' is  the  RMS  velocity  associated  with  scales  of  size  a-jy.  This 
simple  model  is  limited  to  the  early  stages  of  plume  development  where 
a-jy  < Ay.  Once  the  instantaneous  plume  grows  larger  than  the  integral 
scale  the  large  eddies  will  contribute  to  its  growth  and  it  will  cease 
meandering.  In  the  asymptotic  limit  far  downwind  from  the  source, 
a-jy  = ay  so  that  K-jy  = Ky  with 


K.  = A-,  v'A 
iy  7 V 


for  a.  >>  A 

ly  V 


(A22) 


Roth  A5  and  Ay  should  be  constants  of  order  unity.  Using  (A22)  in 
(A20),  and  invoking  the  Lagrangi an-Eul eri an  relation  in  (A6)  leads  to 


da.  f\  . 

^ Lv  for  aiy  » Ay 

dx  Ai  []q 

where  x = U^t  has  been  used  to  convert  from  time  of  travel  to  distance 
coordinates.  By  differentiating  the  exact  asymptotic  solution  (A4)  with 
respect  to  x,  it  is  easy  to  show  that  Ay  = 2 Ai.  With  typical  values 
of  A^  ranging  from  0.35  to  0.80  the  eddy  diffusivity  constant  Ay  is  very 
close  to  unity,  as  expected. 

In  the  early  stages  of  plume  development  meandering  causes  the 
i nstantaneous  and  total  spreads  to  differ.  Using  (A21)  in  (A20)  the  early 
growth  is  given  by 


da . 

]y_ 

dt 


2 Vi'  V 


Expanding  the  derivative  to  2a^^  (da^.^/dt)  gives  the  growth  rate  of  the 
instantaneous  plume  as 


da . 


dt 


V,' 


for  a.y  < 


(A24) 
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with  the  constant  A5  determined  by  the  width  and  shape  of  the  spectral 
bandpass  on  v'  that  produces  vi ' . With  a spectral  bandpass  centered  on  a 
turbulence  scale  of  = a.  , the  effective  velocity  v. ' from  (A15)  can  be 
combined  with  (A24)  to  write,  using  x = U^t 


(A25) 


c V 

This  equation  is  restricted  by  (A24)  to  and  by  (A15)  to  the 

inertial  subrange  of  turbulent  spectrum  where  <<  Because  the 

spectrum  departs  gradually  from  its  inertial  subrange  form,  (A25)  becomes 
progress! vely  less  accurate  as  approaches  A^.  For  homogeneous  turbu- 
lence, (A25)  can  be  directly  integrated  using  the  initial  condition  that 
aiy  - as  at  X = 0 to  obtain 


a.  2/3  . 

ly 


(^)  - 
U A 
C A 


1/3 


a^2/3 


(A26) 


V 

The  constant  A5  may  be  determined  by  comparing  this  result  to  the  result 
of  Smith  (1968)  from  (All)  which  applies  for  immediate  distances  and 
negligible  source  size.  Using  (A6)  to  relate  the  Lagrangian  scale  Ly  to 
the  Eulerian  scale  Ay  and  then  equating  (All)  to  (A26)  yields 


(A27) 


Measured  values  of  Aj  fall  in  the  range  from  0.35  to  0.80,  which  limits 
the  range  of  A5  from  1.4  to  1.9,  which  meets  our  expectation  of  a value 
near  unity.  Using  (A27)  in  (A26)  gives  the  final  form  for  the  instant- 
aneous spread  as 


. a . , 
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1/3 


[( 


3A, 


U 
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-2/3 
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(A28) 
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This  shows  the  expected  a x region  of  accelerated  growth  for  = 0. 

However,  its  most  remarkable  feature  is  that  source  size  effects  do  not 
simply  add  as  a sum  of  variances  for  the  instantaneous  plume.  The  form  of 
(A28)  increases  the  downwind  distance  over  which  source  size  has  a 
significant  effect  on  the  instantaneous  plume  width.  The  source  size  as 
contributes  10%  of  the  instantaneous  plume  width  at  a distance  of 


A u 


JL-  = 15  (-JC) 


(A29) 


s s 

which  for  typical  atmospheric  values  of  v'/Oc  0.10,  Ay  - 100  m and 
Os  - 2 m is  X - 1000  m downwind. 

The  major  deficiency  of  (A31)  is  that  it  assumes  that  plume  spread  is 

always  controlled  by  the  inertial  subrange,  and  continues  indefinitely  at 
3 /2 

an  X growth  rate.  This  form  is  incorrect  far  downwind,  certainly  before 

a.  = A , or  when  a.  = a and  all  meandering  has  ceased. 
iy  V*  iy  y ^ 

However,  in  atmospheric  diffusion  this  limiting  distance  may  be  large, 

because  both  Pasquill  (1974)  p.  228  and  Gifford  (1977)  find  no  evidence  of 

0 5 

cloud  growth  approaching  a « x * asymptote.  This  suggests  that  factors 
such  as  crosswind  shear  and  non-s ta tionari ty  cause  the  apparent  length 
scale  Ay  to  increase  with  plume  travel  time  so  that  meandering  continues 
to  influence  plume  dispersion  at  distances  up  to  several  hundred  kilometers 
in  the  atmosphere. 


Instantaneous  Plume  Spread  Outside  the  Inertial  Subrange 

As  the  instantaneous  plume  width  increases  it  will  be  influenced  by 
large  eddies  at  low  frequencies  outside  the  inertial  subrange,  and  this 
sets  the  limit  of  a-jy  <<  Ay  on  (A29).  To  extend  the  spectral  bandpass 


model  further  downwind  a new  form  for  V)^ ' is  needed.  To  extend  (A15)  to 
the  large  eddy  end  of  the  spectrum,  we  assume 
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V ( 


1/3 


Z.  + A ' 
k V 


(A30) 


which  produces  the  correct  value  of  Vj^'  ->  v'  as  Using  the 

definition  of  the  spectral  density  E2(k)  in  (A17)  and  scale  J?.j^  in 
(A16),  the  above  equation  leads  to  a spectrum  of  the  form 


2 A.A  v'^ 

E^Ck)  = i-J! — (A31) 

3{1  + 

which  has  the  correct  integral  value 

00 

v'2  E / E2(k)  dk  (A32) 

0 

for  any  value  of  A5.  The  value  of  A5  is  determined  by  the  zero 
frequency  intercept  which  must  be  equal  to 

E2(0)=iv'^A  (A33) 

which  sets  the  constant  at  A^  = 3/ti.  With  this  value,  the  spectrum  pro- 
duces an  inertial  range  spectrum  at  high  wavenumber  of 

E^  (k)  = 0.69  v'^  a'^^^  k‘^^^  (A34) 

The  constant  0.69  is  in  good  agreement  with  measurements  in  the  inertial 
subrange  at  high  Reynolds  number,  see  Hinze  (1975)  p.  255.  This  result 
leaves  us  with  some  confidence  in  (A30)  as  a reasonable  representati on  of 
the  spectral  bandpass  function.  Using  (A30)  in  (A24)  with  a spectral 


bandpass  centered  at  = o]y 
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a.  + A 1/3  , 

da.  = U^]  dx 
^ ly  6^u  ^ 

c 


a . 

ly 


To  produce  a suitable  closed  form  integral  we  expand  the  left  side  of 

a . + A 1/3  A 1/3  a . 1/3 


(• 
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a. 

ly 


a.  A 

iy 


and  approximate  the  second  term  to  write 


a.  .A  A 1/3  2/3 


aiy 


aiy 


4 ^A, 


This  approximate  form  is  within  10%  of  (A36)  in  the  range  0 < 
Integrating  (A35)  using  (A37),  from  a^  at  x = 0 to  at  x,  with  A^ 
evaluated  from  (A27) 

a 2/3  a 4/3  1/3 

( JZ)  + 1 ( JZ)  = (_^)  (Z_)  Z.  + $ 

Z ^ 8 Z ^ ^3A/  ^ A s 

V V lev 


wi  th 


a^  2/3  , a^  4/3 

$ = (_i)  +1  (_i) 

s ^A  ^ 8 ''A 

V V 


Equation  (A38)  is  a quadratic  with  the  positive  root 


a.  2/3  1/3  , 1/2 

[_iy]  =[16+8((J_]  + -4 

A 3A,  IJ  A 

V lev 


(A35) 

(A35) 

(A36) 

(A37) 
< 30. 

(A38) 

(A39) 

(A40) 


which  will  be  valid  for  a-jy  < Ay,  a limit  imposed  by  (A24). 


APPENDIX  B 


ALONG-WIND  DEVELOPMENT  OF  CONCENTRATION  FLUCTUATIONS 


In  Wilson,  Robins  and  Fackrell  (1982)  a model  for  along  wind  changes  in 
fluctuation  intensity  was  developed,  based  on  a balance  between  along-wind 

advection  and  dissipation  of  c'^.  The  effects  of  source  size  were  dealt  with 
empirically  through  a virtual  origin  that  depended  strongly  on  source 
di ameter. 

The  major  deficiency  of  this  model  was  its  tendency  to  underpredict  the 
intensity  of  fluctuations  close  to  an  elevated  source.  This  underprediction 
is  caused  by  the  assusmption  that  the  plume  spread  parameter 
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is  proportional  to  x * . This  will  only  be  true  very  far  downwind,  and  near 
the  source  the  relation  is  close  to  X « x.  In  the  following  analysis  we  will 
correct  this  deficiency,  and  at  the  same  time  make  use  of  a new  way  of 
scaling  dissipation  that  accounts  explicitly  for  source  size  effects. 

Using  the  experimental  observation  that  dissipation  and  advection  are  in 
approximate  balance  at  each  point  in  the  plume  we  have 


U .Ac„  = 
^ dx 


”£ 

c 


(B2) 


where  U^-  is  the  convection  velocity  of  the  plume,  and  eq  is  the 
dissipation  of  concentration  fluctuation  variance,  which  we  will  model 
through  the  dissipation  time  scale  Tg- 
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Assuming  that  the  variance  profiles  have  the  self-preserving  form 


c'2  = c g(i-,  2-)  (R4) 

a a 

y z 

we  can  combine  (B2),  (B3)  and  (B4)  and  integrate  across  the  plume  to  obtain 

(B5) 
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^ dx 
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Because  source  height  is  not  a parameter  this  result  must  be  limited  to  plumes 
with  h » 02  or,  to  ground  level  sources  with  h = 0 for  which  self- 
preserving similarity  solutions  also  apply. 

The  integrated  balance  equation  (B5)  was  developed  from  an  advection- 

dissipation  balance.  Sykes,  Lewellen  and  Parker  (1984)  developed  an  exact 

~ —2  — 2 

relation  by  writing  the  balance  equation  for  the  square  c = c + c'  . The 
mean  square  contains  no  production  terms,  and  the  diffusion  terms  vanish  when 
integrated  over  the  entire  plume.  For  the  special  case  of  Gaussian  profiles 
for  both  "c  and  c'^  they  obtain,  in  the  appendix  to  their  paper 


c2 

m- 


dx 


f , I 

(c  + :r- 


(B6) 


where  the  c^  is  the  spatial  average  of  the  mean  concentration  over  the  plume 

_ 2 

cross  section  (i.e.  in  y and  z).  The  c^  /2  results  from  squaring  the  Gaussiar 

— — 2 
for  c before  integration.  Their  results  show  that  the  c^  /2  term  is  only 

important  very  close  to  the  source,  and  at  all  distances  of  practical  interest 

(B6)  produces  the  same  results  as  (B5). 
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Modelling  the  Dissipation  Time  Scale 

One  important  contribution  of  Sykes,  Lewellen  and  Parker  (1984)  was  to 
suggest  a physically  realistic  model  for  the  dissipation  time  scale.  They 
suggest  that  should  be  related  to  the  properties  that  influence  the  dev- 
elopment of  the  instantaneous  plume,  rather  than  the  total  plume  spread  which 
includes  meandering.  The  reasoning  behind  this  is  that  because  meandering 
acts  only  to  push  the  entire  plume  around  it  does  not  contribute  to  dissipa- 
tion, which  is  dominated  by  small  scale  motions.  Using  the  instantaneous 
spread  a-j  and  the  RMS  turbulence  velocity  scale  u^-j  of  eddies  of 
comparable  size  to  the  plume  we  define 


(B7) 


Then,  using  the  spectral  scaling  developed  in  Appendix  A,  the  velocity  scale 
of  eddies  of  size  a-j  must  be  related  to  the  total  turbulence  scale  by 
(A15).  For  growth  in  the  inertial  subrange  of  the  turbulence  spectrum 


so  that  (B6)  becomes 
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Following  the  turbulent  diffusivity  model  for  instantaneous  spread  from 
Appendix  A,  we  can  express  plume  spread  a-j  in  a form  similar  to  (A24) 


da . 
1 

dx 


(BIO) 
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As  shown  in  Appendix  A this  equation  is  valid  only  for  a-j  < A.  Using  (B8), 
integrate  (RIO)  to  obtain  equations  similar  to  (A26)  and  (A28)  in  Appendix 
A to  obtain 

2/3  2/3 

(^5  (BID 

A ^ U A A 

C 

as  suggested  by  Sykes,  Lewellen  and  Parker  (1984).  Then,  assume  the  total 
spread  a,  near  the  source  is 

ik  = ^ (B12) 

dx  U 

c 

so  that 


(B13) 


Here  we  have  made  the  assumption  that  the  initial  source  size  has  a 
negligible  effect  on  the  total  spread,  and  affects  only  the  instantaneous 
plume.  Combining  (B9),  (Bll)  and  (B13)  we  find  the  dissipation  time  scale  is 
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Note  that  the  influence  of  source  size  is  accounted  for  explicitly  by  the 
virtual  origin  term  a^/A, 

At  this  point  we  must  decide  on  an  appropriate  definition  for  the  plume 
spread  a,  and  the  turbulence  scales  Ug  and  A.  The  process  of  dissipation 
of  c'^  is  three  dimensional  in  nature  and  is  isotropic  in  the  smallest 
scales.  With  this  in  mind,  the  form  of  the  integrated  balance  equation  (B5) 
suggests  the  form 
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With  nothing  to  guide  us  in  our  definition  of  the  turbulence  length  scale  A, 
we  arbitrarily  set 


A 


A + A 
V w 


(B16) 


By  transforming  the  variance  balance  equation  (B5)  from  x to  a coordinates, 
we  can  avoid  the  problem  of  defining  the  turbulence  velocity  scale  U5. 

Using  (B12) 


A-  = H3.A- 

dx  U da 
c 

Then  with  (B17),  (B15)  and  (B14),  we  can  transform  (B5)  to 
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The  transformed  balance  equation  may  be  integrated  to  obtain 

c S = 


S + S ] 3 

0^ 


(B19) 


This  can  then  be  normalized  by  the  mean  concentrati on  Cq  on  the  plume 
centerline,  using  a plume  mass  balance  to  show  that 
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which  is  constant  with  downwind  distance.  Defining  another  normalized 
distance  as 

. ^n.5 


X = 


a o 

y z 
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(B21) 
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where  H is  a convenient  constant  length  scale.  In  homogeneous  turbulence 
H E A,  while  in  shear  f 1 ow  H may  be  taken  as  the  shear  layer  thickness. 

Using  this  distance  parameter  the  fluctuation  intensity  i = c/T^  may  be  found 
from  the  ratio  of  (B19)  to  (B20) 


i = B,  {-]  


(B22) 


where  a 
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B2/2,  and  Bg  = (B^/B^)  The  virtual  origin  is 
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By  differentiating  (B22)  it  is  easy  to  show  that  the  intensity  i has  a maximum 

value  at  X = X /(a-1). 

max  s ' 

The  initial  increase  of  i near  the  source  is  due  to  the  different  decay 

" _ 2 
rates  of  c and  c^ , and  not  to  production  of  c'  as  might  seem  to  be  the 

2 

obvious  cause.  In  fact,  production  of  c'  is  neglected  in  the  advection- 

dissipation  balance.  Close  to  the  source  where  S <<  the  RMS  fluctuation  in 

(B19)  decays  like  c ~ S”\  while  far  downwind  where  S » S^,  c ~ S 

— -2 

Because  the  mean  in  (B20)  ~ S at  all  locations,  the  ratio,  c/c^  will  at 

first  increase  as  i ~ S and  finally  decrease  as  i ~ S^”^,  even  though  both  c 

and  c"  decrease  continuously. 

0 
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Equation  (B22)  for  the  along-wind  development  of  intensity  i has  an  im- 
plicit limitation  to  homogeneous  turbulence,  where  A is  constant  and  there 
are  no  ground  surface  dissipation  effects.  In  the  next  section  we  will  relax 
the  condition  of  homogeneous  turbulence  and  use  (B22)  and  (B23)  for  both 
elevated  and  ground  level  sources,  where  the  effective  scale  A must  be  re- 
garded as  a distance  dependent  average  over  the  plume  travel. 


Effective  Turbulence  Scale  A in  Shear  Flows 

If  we  consider  the  limiting  cases  of  a very  high  source,  and  a ground 
level  source  we  find: 

• For  a high  source  where  h » the  turbulence  is  approximate- 
ly homogeneous,  and  the  length  scale  A should  be  evaluated  at 

z ^ = h. 
ref 

• For  a ground  level  source  the  plume  encounters  scales  which 
increase  with  downwind  distance  as  the  centroid  of  the  concen- 
tration profile  increases  its  distance  above  the  surface. 

Because  Lagrangian  similarity  theory  predicts  that  both  the 
centroid  and  the  plume  spread  increase  linearly  with  time 
of  travel , it  is  reasonable  to  evaluate  the  effective  scale  A 
at  = Bga,  where  Bg  is  a constant. 

To  interpolate  between  these  limiting  cases  we  define  a reference  height 

2 2 

MRfiO)  ] (B24) 


Finally,  assuming  a simple  linear  relation  between  the  effective  scale  and 
the  reference  height 
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Bpj  z . 
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Intensity  of  Conditionally  Averaged  Fluctuations 

The  conditionally  averaged  intensity  i^  has  the  intervals  of  zero 
concentration  removed,  and  at  the  same  time  has  no  ground  surface  effects 
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acting  to  dissipate  fluctuations.  To  derive  a conditional  "plume"  intensity 
that  is  realistic  in  inhomogeneous  turbulence,  define  i^  as  the  intensity  in 
the  distance  dependent  turbulence  experienced  by  a ground  level  source,  but 
with  no  surface  dissipation  effects.  Although  this  hypothetical  case  is 
difficult  to  simulate  experimentally,  it  may  be  derived  analytically  by 
setting  h = 0 in  (B24)  to  find  A (R25),  and  using  this  effective  scale  in 
(B22)  to  obtain 
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where  from  (R24)  and  (B25)  and  the  definition  of  \ in  (B21) 


and  from  (B23) 
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This  simple  method  for  deducing  the  conditionally  averaged  intensity  is  based 
on  the  observation  that  there  is  no  plume  meandering  or  small  scale  intermit- 
tency  on  the  plume  axis  for  a ground  level  source.  Under  these  conditions, 
where  diffusion  occurs  in  a field  of  very  small  scale  turbulence,  the  condi- 
tionally averaged  and  total  plume  statistics  are  identical  because  the  inter- 
mi  ttency  factor  is  unity. 


APPENDIX  C 


ENTRAINMENT  AND  DILUTION  RY  SELF-GENERATED  TURBULENCE 

The  theory  necessary  for  estimating  the  dilution  in  the  buoyant  jet  by 
turbulence  generated  within  the  plume  is  developed  in  this  Appendix.  This 
sel f-generated  turbulence  is  caused  by  the  momentum  and  buoyancy  fluxes  from  the 
source,  and  causes  entrainment  of  ambient  air  into  the  plume.  Close  to  the 
source,  this  self-generated  entrainment  is  the  most  important  mechanism  for 
dilution.  As  the  plume  rises  and  is  carried  downwind  the  momentum  jet 
turbulence  rapidly  decays  and  buoyancy  generated  turbulence  becomes  dominant. 
Finally,  atmospheric  turbulence  becomes  the  dominant  factor  in  plume  dilution  as 
the  decaying  buoyancy  generated  turbulence  is  overpowered  by  the  constant 
atmospheric  turbulence. 

On  the  following  sections  the  initial  dilution  in  the  initial  momentum  jet 
will  be  estimated.  Then,  the  dilution  due  to  combined  buoyancy  and  momentum  in 
the  bent  over  plume  will  be  computed  in  terms  of  the  plume  rise. 

Dilution  in  the  Source  Jet 

To  simplify  the  analysis,  assume  that  the  jet  has  travelled  far  enough  to 
entrain  sufficient  air  that  its  density  is  close  to  the  ambient  air  density 
p0.  This  distance  will  depend  on  the  composition  of  the  pipe  gas,  but  will  be 
at  least  as  far  as  the  distance  to  the  point  where  compressibility  effects  are 
negligible.  Wilson  (1981  b)  pp.  36-37  showed  that  this  is  about  x > 30  Dg, 
where  Dg  is  the  source  exit  diameter. 

The  experimental  results  of  Ricou  and  Spalding  (1961)  showed  that  the 
process  of  air  entrainment  is  kinematic  rather  than  dynamic,  so  that  the 
entrainment  rate  depends  only  on  jet  velocity  and  not  on  its  density.  The 
entrainment  velocity  is  defined  by 
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V = a W 
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with  an  entrainment  constant  a that  is  independent  of  densities  pa  and  pe. 
Because  the  entrained  mass  is  at  density  pg,  the  mass  added  to  the  jet  in 
distance  dz  is 

dm  = p^v^  (2ixRdz)  (C2 

Using  (Cl)  this  can  be  written 

^ = 2n  a p R W (C3 

dz  ® 

Because  both  jet  velocity  and  radius  vary  with  distance,  a momentum  balance  mu 
be  used  to  relate  R and  W.  Using  the  approximation  that  the  local  jet  density 
is  equal  to  the  ambient  air  density  p^  allows  us  to  write  the  momentum  balance 
in  the  constant  pressure  jet,  which  begins  after  its  initial  expansion  from  "e 
to  "s",  (the  "source"  condition) 


so  that 


p 11  R^w2  = p u r2w2 
s s 
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Because  there  is  no  crater  drag  interaction  here,  the  "source"  condition  is 
equivalent  to  the  "expanded  jet"  condition  denoted  by  a subscript  "j"  in  Wilsc 
(1981  b).  Using  (C5)  in  (C3), 


M = 2_a  r^i 
dz  R^  p^ 
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where  m = m = p'nR  W is  the  mass  release  rate.  Integrating  this  between 
e s *^s  s s ^ ^ 

the  origin  at  "s"  and  a height  z yields  the  mass  dilution  ratio 
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The  area  ratio  A /A  = R /R  has  been  introduced  to  allow  distance  z to  be 
e s e s 


normalized  with  the  actual  pipe  exit  radius  R^  rather  than  the  hypothetical 
source  radius  R^. 

Dilution  in  a Rising  Bent-Over  Plume 

Briggs  (1975)  shows  that  the  transition  from  a vertical  jet  to  a bent  over 
plume  occurs  very  rapidly,  and  that  most  of  the  trajectory  can  be  approximated 
as  fully  bent  over,  with  the  plume  axis  nearly  parallel  to  the  wind.  As  the 
plume  bends  over  a pair  of  counter  rotating  vortices  develop  which  enhance 
entrainment  into  the  underside  of  the  plume.  At  the  same  time,  the  vertical 
velocity  W moves  the  bent  over  plume  perpendicul ar  to  its  axis,  causing  more 
entrainment  than  in  a vertical  jet.  To  account  for  this  we  define  an  enhanced 
entrainment  coefficient  p. 


Spillane  (1983)  reviews  several  other  entrainment  models  that  can  be  applied  to 
the  transitional  bending-over  phase.  Typically  p = 0.4  to  0.8,  while  the 
vertical  jet  entrainment  coefficient  is  only  a = 0.08.  By  the  time  the  jet  is 
bent  over,  it  has  entrained  enough  ambient  air  that  its  density  is  close  to 
Pa.  Because  the  convection  velocity  is  a constant  U^,  rather  than  the  jet 
velocity  W that  varies  within  z,  a mass  balance  yields 


(C8) 


dz 


(C9) 


Defining  an  effective  initial  radius  Rq  that  contains  the  mass  emission  rate 
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m^,  see  Figure  4-1,  and  integrating  (C9)  from  z = 0 to  the  plume  rise  height 


z = Ah 


R = R + pAh 
0 


(CIO 


Because  the  plume  has  a density  pg  and  is  carried  along  at  speed  the 

initial  radius  R must  be  related  to  m by 
0 e 


m = p u R U. 
e ^a  0 h 


(Cll 


The  total  mass  flux  through  a plume  cross  section  is 


m = p u R IJ. 
"^a  h 


(C12 


Taking  the  ratio  of  (C12)  to  (Cll)  and  using  (CIO)  the  mass  dilution  ratio  is 
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^ = (1  + g...  Ah) 
m R 
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At  the  pipe  exit  the  mass  flow  is 


m = p R W 
e ^e  e 


(C14 


Equating  (Cll)  and  (C14)  allows  the  effective  radius  Rq  at  the  origin  to  be 
expressed  in  terms  of  the  pipe  exit  radius  Rq 
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Using  this  in  (C13)  the  mass  flow  in  the  bent  over  plume  is 
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Compare  this  to  the  mass  dilution  for  a vertical  jet  given  in  (C7).  For  a jet 
the  dilution  is  linear  in  rise  distance,  while  for  a bent  over  plume  it  varies 
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as  the  square  of  the  plume  rise.  The  more  rapid  dilution  of  a bent  over  plume 
is  due  to  the  constant  convection  velocity  U^,  in  contrast  to  the  convection 
velocity  W which  decreases  with  distance  in  a vertical  momentum  jet. 

High  Pressure  Pipeline  Ruptures 

The  mass  dilution  ratios  for  a vertical  jet  and  a bent  over  plume  in  (C7) 
and  (C16)  apply  to  any  fluid,  as  long  as  the  jet  is  sufficiently  diluted  to  have 
a density  near  ambient  p^.  In  this  section  the  flow  character! sti cs  and  fluid 
properties  of  a high  pressure  choked  (i.e.  sonic  velocity)  gas  release  will  be 
used  to  express  the  dilution  in  terms  of  the  readily  available  pressure  Pg  and 
stagnation  temperature  T^  at  the  pipe  exit.  From  the  ideal  gas  law  the 
density  ratio  in  (C7)  is 


where  A is  the  molecular  weight  and  T the  absolute  temperature.  The  molecular 
weight  remains  the  same  during  the  expansion  because  no  air  is  entrained.  With 
Ai  = Ae  = A^  we  will  use  A^  throughout  as  the  pipe  gas  molecular  weight.  For 
adiabatic  flow  the  jet  temperature  T5  is 


where  is  the  Mach  number  in  the  expanded  jet.  Wilson  (1981b)  p.  19  shows 
that  for  high  pressure  at  the  pipe  exit 


(CI8) 
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as  long  as  Pe/Pa  > 10.  Using  (C19)  in  (C18) 
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(C20) 


Wilson  (1981b)  p.  19  also  showed  that  even  at  moderate  exit  pressure  ratios. 


A 1 .P  > 
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By  using  (C17),  (C20)  and  (C21)  in  (C7)  the  mass  dilution  ratio  for  a vertical 
jet  becomes 


where 


a-  = 1 + 2 a a,(l) 

• 1^  L ^ 
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The  approximation  in  (C22)  is  limited  to  high  exit  pressures  Pe/^a  ^ 
to  distances  far  enough  from  the  exit  that  the  jet  density  is  approximately 
equal  to  that  of  the  ambient  air. 

Next,  consider  the  equivalent  process  for  a bent-over  plume  whose  mass  f 1 o 
is  given  by  (C16).  To  evaluate  the  density  ratio  and  exit  velocity  in  (C16)  th 
choked  flow  condition  of  = Sg  is  used,  where  Se  is  the  speed  of  sound 
at  the  pipe  exit  temperature 
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The  first  bracketed  term  is  the  speed  of  sound  s\  in  the  pipe  gas  at  its 
stagnation  temperature  The  temperature  ratio  in  the  second  term  is  given 
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where  the  Mach  number  Mg  = 1.0  in  the  choked  exit  flow,  so 

T. 


1 = k+1 

T 2 
e 


Using  this  in  (C24) 
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From  the  ideal  gas  law  the  density  ratio  in  (C16)  is: 
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with  (C26)  and  (C27)  this  allows  us  to  write  (C16) 
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where  the  constant  a 


^ is  defined  by  (C23) 
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Dilution  Factors  for  Mass  and  Volume  Concentration 

The  concentration  of  a toxic  or  flammable  component  in  a mixture  is  often 
expressed  in  terms  of  its  volume  fraction,  c.  By  assuming  that  this  component 
and  the  entire  mixture  both  behave  as  ideal  gases  we  can  develop  the  relation- 
ship of  volume  fraction,  c,  to  the  mass  fraction  f.  It  is  easy  to  show  that  the 


dilution  factor  for  mass  fraction  is  just  the  ratio  of  mass  flow  rates.  That 
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where  f is  the  mass  fraction  of  the  component  at  distance  z from  the  exit,  and| 
f^  is  the  mass  fraction  of  this  component  in  the  pipeline  gas.  This  mass 
fraction  dilution  factor  may  be  directly  determined  from  (C22)  for  a jet,  or 
from  (C29)  for  a bent  over  plume. 

The  relationship  of  mass  and  volume  fraction  may  be  found  from  the  ideal 
gas  law.  Using  the  law  of  partial  volumes,  the  volume  fraction  c and  mass 
fraction  f of  a specific  component  gas  are  related  by 


c = f .jsnx 

A 


(c: 


where  A is  the  molecular  weight  of  the  component  and  Af^-jx  molecular  weigh' 
of  the  mixture. 

During  the  dilution  process  air  is  entrained  by  the  jet,  changing  the 
molecular  weight  of  the  mixture  from  A^ix  = Ag  at  the  pipe  exit  point  to 
^mix  ^ ^a  molecular  weight  of  ambient  air.  Taking  the  ratio  of  the 
volume  fraction  c^  in  the  pipeline  gas  to  the  jet  concentrati on , c,  at  a 
point  far  enough  from  the  exit  to  have  a density  close  to  ambient  air,  we  finj 
the  volume  dilution  factor,  c^/c 
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Using  (C30) 


c A • 
e e r m 


c A 


(C33) 


ni 

3 e 


Combine  (C29)  and  (C33),  with  (C23)  to  evaluate  the  parameter  a^,  and  we  find 
for  a bent  over  plume 


A f-T,  P 

_^)  + P fJi  _i 

A TP 

a U a e 


0.5  0.5 


2 

k+1 


IJ 


and  for  a jet,  using  (C22),  (C23)  and  (C33)  the  dilution  factor  is 


_ + 2 a 
A 

a 


0.5  n 
z 

F 


(C35) 


It  is  interesting  to  note  that  for  a horizontal  half-conical  jet  bounded  by  the 
ground,  the  change  in  geometry  will  affect  both  the  entrainment  and  momentum 
conservation  in  (C2)  and  (C4)  in  the  same  way.  Because  of  this,  the  dilution  in 
a "half  jet"  will  also  be  given  by  (C35). 


Transition  From  Vertical  Jet  to  Bent-Over  Plume 

As  the  wind  bends  the  vertical  jet  over,  its  mass  flow  rate  dilution  will 
gradually  change  from  (C7)  to  (C16).  We  will  assume  that  the  point  of 
transition  where  (C16)  can  be  applied  occurs  when  the  plume  is  bent  over  at  an 
angle  of  45°.  At  this  point,  the  plume  trajectory  will  have  a slope  of  unity. 

By  simple  kinematics  the  vertical  velocity  W of  the  jet  must  equal  the  windspeed 
U^'  that  is  bending  the  plume  over.  Setting  = W in  (C5) 


W 

U 


1 

h 


(36) 


8 


Wilson  (1981b)  suggests  that  density  differences  will  have  little  effect  on  je 
spreading  rate,  so  that 


(C3 


To  simplify  the  transition  condition,  assume  that  the  virtual  origin  of  the  je 
is  at  the  pipe  exit.  Integrating  (C37)  using  this  assumption  of  R = 0 at  z = 


R = 2 a z 

so  that  transition  to  bent-over  plume  behavior  in  (C36)  occurs  when 


(C3) 


R 


'“  Pa  ^'h 
transi ti on 


(c: 


This  is  only  an  approximate  estimate  because  the  point  of 
arbitrarily.  Also  the  entrainment  "constant"  a increases 
bending  over,  making  a in  (C39)  difficult  to  define. 


45°  slope  was  choser 
during  the  process  c 


APPENDIX  D 


SUM  OF  CURES  APPROXIMATION  TO  FINAL  RISE 


In  this  appendix  Briggs  (1975)  plume  rise  model  will  be  reviewed  by 
presenting  the  theoretical  equations  for  plume  rise  due  to  the  combined 
effects  of  stack  gas  buoyancy  and  vertical  jet  momentum.  Briggs'  turbulent 
dissipation  criterion  for  determining  the  point  of  final  rise  depends  on  an 
estimate  of  the  rate  of  dissipation  F within  the  plume  due  to  its  self- 
generated turbulence.  His  estimate  will  be  shown  to  be  consistent  with  the 
local  scaling  model  developed  in  Appendix  A.  The  times  of  travel  to  the 
point  of  final  rise  of  pure  momentum  jets  and  pure  buoyant  plumes,  will  be 
compared  to  the  travel  time  of  a combined  buoyant  jet  to  its  final  rise. 

The  combined  momentum  and  buoyant  rise  of  a plume  depends  strongly  on 
the  rate  of  entrainment  of  ambient  air  into  the  plume.  This  is  characteri zed 
by  an  entrainment  velocity  v^,  which  for  a bent-over  plume  with  its  axis 
close  to  horizontal  is  related  to  the  vertical  local  rise  velocity  W 

v^  = PW  (Dl) 


as  shown  schematically  in  Figure  4-1.  Briggs  (1975)  solved  the  mass, 
momentum  and  thermal  energy  conservation  equations  for  the  case  where  the 
entrainment  "constant"  p is  independent  of  x,  (but  may  still  be  a function  of 
the  same  source  velocity  ratio  W^/U^).  In  this  case  the  plume  rise  Ah 
reduces  to  the  simple  closed  form 


Ah 


^ t + — — 


(D2) 
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in  terms  of  the  plume  travel  time,  t.  Equation  (D2)  also  requires  a uniform 
wind  speed,  constant  with  z throughout  the  plume  layer  from  z = h^  to 
z = h^  + Ah.  This  constant  wind  speed  allows  the  travel  time  t to  be 
related  to  downwind  distance  by  x = U^t.  The  key  parameters  in  (D2)  are  the 
momentum  and  buoyancy  fluxes,  and  F,  which  are  defined  as 


(n3) 


and 


F = g 


(— 


.1  W R' 

} c < 


(D4) 


Final  Rise  Using  Briggs  Criteria 

Briggs  (1975)  assumes  that  final  rise  occurs  when  the  average 
dissipation  rate  "e  of  self-generated  turbulence  energy  within  the  plume  is 
equal  to  the  local  atmospheric  turbulence  dissipation  rate  F.  When  these  tw 
rates  are  equal  Briggs  hypothesizes  that  the  atmospheric  turbulence  is  able 
to  "invade"  the  plume  and  decorrelate  its  mixing  processes,  causing  it  to 
"break  up".  To  make  use  of  this  condition  to  predict  final  rise,  we  must 
first  express  the  internal  dissipation  rate  in  terms  of  plume  rise.  Using 
large  eddy  scaling  of  turbulent  dissipation  (see  (A13)  in  Appendix  A) 


(05 


t 


where  w^  is  the  root  mean  square  turbulence  velocity  within  the  plume,  and  l| 
is  the  scale  of  large  self-generated  eddies  within  the  plume.  The  constant 
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- 0.8  from  experimental  measurements  of  a wide  range  of  turbulent  flows. 

If  we  make  the  rather  crude  estimate  that 

w‘  - W (D9) 

P 

1 = R (D7) 

P 

where  R = pAh  is  the  plume  radius.  We  obtain,  for  = 0.8  and  p = 0.6 

e = 1.33  _ (08) 

Ah 

Rather  than  constrain  ourselves  with  this  estimate,  we  will  follow  Briggs* 
convention  and  leave  the  constant  to  be  determined  later  as 


(09) 


Briggs  (1975)  p.  75  estimates  that  ti  = 1.5,  which  is  quite  close  to  the 
independently  computed  value  of  ti  = 1.33  in  (08).  The  time  of  travel  until 
final  rise  is  found  by  observing  that 

W = (DIO) 

dt 


and  using  the  combined  rise  equation  (02)  to  evaluate  both  W and  Ah  in 
(010).  Setting  internal  and  external  dissipation  rates  equal  at  the  point 
where  final  rise  occurs,  T = e we  find  after  some  manipulation 


+ 1 

.2  U. 


9/7 


2 IJ. 


3/7 


(Oil) 


where  t^  = x^/U^  is  the  plume  travel  time  at  the  point  of  final  rise. 


90 


Final  Rise  of  Pure  Buoyant  Plumes  and  Momentum  Jets 

For  a pure  momentum  jet  the  buoyancy  f 1 ux  F = 0,  and  (Dll)  gives  the 
time  to  final  rise  as 


(D12) 


and  a final  rise  height  from  (D2)  with  F = 0 of 


3 '"m  ; 
>if 


1/3 


(D13) 


The  superscript  " ~ " indicates  a pure  momentum  jet  or  buoyant  plume.  For  a 
pure  buoyant  plume  the  momentum  flux  F^^  = 0 and  from  (Dll)  and  (02)  the  time 
to  final  rise  tpf  and  the  rise  height  Ahef  are 


and 


Inserting  (D12)  in  (D13)  and  (D14)  in  (D15)  leads  to  the  form  given  by  Briggs 
(1975) 


(D16 


and  for  the  buoyant  plume 


Ah 


Rf 


L3p^  hJ  ^ 


{D17 


91 

It  is  clear  from  these  expressions  that  variation  with  height  of  the 
atmospheric  turbulence  dissipation,  e,  must  be  known  in  order  to  determine 
the  final  rise. 


Variation  of  Atmospheric  Dissipation  with  Height 

Most  stack  designs  are  based  on  dispersion  in  neutral  atmospheric 
stability,  where  mechanically  generated  turbulence  from  wind  shear  is  the 
only  mechanism  for  dispersion.  Neutral  stability  produces  a log-law  velocity 
profile  in  the  surface  layer  so  that 


U = ^ In  (^) 
1C  z 

0 


(ni8) 


where  k = 0.4  is  von  Karman's  constant.  It  can  be  shown  that  turbulent 

dissipation  in  this  log-law  surface  layer  is 

3 

e = (D19) 

K Z 

At  plume  height,  where  z = (h^  + Ah)  this  becomes 


e 


K(h^  + Ah) 


(D20) 


This  variation  of  e with  height  leads  to  closed  form  solutions  for  plume  rise 
for  the  limiting  cases  of  tall  stacks  and  for  ground  level  releases.  For 
tall  stacks  in  high  winds,  where  h^  » Ah  we  find 


3 

e = = CONSTANT 

< h 

s 


(D21) 
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A constant  value  of  e with  height  is  also  typical  of  unstable  atmospheric 
mixing  where  convective  turbulence  is  dominant,  at  times  of  strong  solar 
heating  of  the  ground  combined  with  light  winds.  The  approximation  for  e in 
(D21)  does  not  apply,  but  the  general  results  for  constant  e in  the  following 
analysi s wi 1 1 be  val id. 

Another  case  of  practical  interest  is  the  release  from  a ground  level 
source  into  mechanical  turbulence.  In  this  case  the  stack  height  = 0 and 

turbulence  dissipation  varies  directly  with  plume  rise  height  Ah, 

3 

e = — — (D22) 

K Ah 

In  the  following  sections  combined  rise  equations  for  these  limiting  cases  of 
£ = constant  and  e « Ah“^  will  be  developed. 

Normalizing  Plume  Rise  Variables  to  Obtain  Parametric  Solutions 

When  the  stack  height  h^  is  much  larger  than  the  plume  rise,  or  in  the 
convective  turbulence,  the  atmospheric  turbulence  dissipation  e is  almost 
constant  through  the  plume  rise  layer.  In  the  first  case  considered  below  we 
will  assume  that  e is  a constant,  and  use  the  ratios  of  plume  rise  so  that  e, 
the  entrainment  constant  p,  and  other  difficult  to  determine  parameters 
conveniently  cancel. 

The  ultimate  objective  of  this  analysis  is  to  determine  the  ratio  of  the 
actual  combined  rise  of  a plume  with  both  momentum  and  buoyancy  to  the 
approximate  rise  found  by  considering  buoyancy  and  momentum  rise  separately 
and  then  adding  them  as  a sum  of  cubes.  The  only  difficulty  is  that  an 
explicit  solution  can  not  be  found  for  this  ratio.  Instead,  several 
equations  are  required,  connected  through  a dimensionless  momentum  flux  to 
buoyancy  flux  ratio  which  is  defined  as 
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^ = 


2 F t, 


Mf 


(023) 


where  F and  Fj^  are  the  buoyancy  and  momentum  fluxes  and  t^^^^  from  (D12)  is  the 
time  to  final  rise  of  a pure  momentum  jet  with  the  same  momentum  flux  and 
entrainment  constant  as  the  buoyant-momentum  plume.  For  the  special  case 
where  the  atmospheric  turbulence  dissipation  e is  a constant  with  height,  the 
mysterious  collection  of  terms  used  to  define  ^ can  be  directly  related  to 
the  pure  buoyant  and  momentum  final  rises  in  (D16)  and  (017)  by 


5 = for  e s CONSTANT  (D24) 

L*^Mf  J 

The  time  of  travel  t^  to  the  point  of  final  rise  for  the  combined  plume 
can  be  expressed  in  terms  of  ^ and  tne  final  rise  time  t^^^  of  a pure  momentum 
jet  by  using  (012)  and  (023)  in  (Oil)  to  obtain 

(1  + L r,)5/7 

r,  = 2 — i for  j.  = CONSTANT  (025) 

where  the  time  ratio  r^  is  defined  as 

t. 

r,  = for  e = constant  (D26) 

^ t 
^Mf 

Note  that  (025)  is  an  implicit  equation  that  must  be  solved  by  iteration  for 
each  specified  value  of  the  rise  ratio 


Similar  equations  can  be  developed  for  the  case  of  a ground  level  re- 
lease into  mechanical  turbulence  where  the  atmospheric  turbulence  dissipation 
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e is  a function  of  height,  given  by  (D22).  Using  (D22)  in  (D12)  to  get  t^^ 
and  then  using  this  in  (023), 

Ahf5^  -1 

^ for  e « (Ah)  (027) 


For  this  case,  the  parametric  variable  ^ is  simply  the  ratio  of  buoyancy  to 
momentum  final  rise.  Then,  using  (023)  in  (Oil)  and  combining  (022)  and 

(ni2). 


'^2  " 


(1  ^Ir^) 


for 


(Ah) 


-1 


where  the  time  ratio  r^  is  defined  in  the  same  way  as  r. 


(028) 


-1 

for  e « (Ah) 


(029) 


Now,  with  the  parametric  variables  r^  and  r^  defined,  we  can  derive 
equations  for  the  ratio  of  combined  final  rise  to  its  sum  of  cubes  approxima- 
tion. 

Ratio  of  Combined  Final  Rise  to  "Sum  of  Cubes  of  Final  Rises*'  Approximation 
The  transitional  plume  rise  equation  (02)  clearly  shows  that  momentum 
and  buoyancy  effects  add  according  to  a sum  of  cubes  law.  This  is  easy  to 
see  when  (02)  is  written  at  the  point  of  final  rise 


Ahf 


(030) 


where  the  buoyancy  and  momentum  rise  components  are,  at  final  rise  where  the 
travel  time  t = t^ 
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(031) 


1/3 


Ah, 


Bf 


(032) 


with  these  definitions  this  sum  of  cubes  law  is  an  exact  expression. 

A simpler  approximate  form  for  the  final  rise  results  when  we  apply  the 


components  as  if  they  were  two  separate  plumes.  These  separate  plume  rises 
have  been  designated  by  a superscript  " ~ " in  (016)  and  (017)  to  distinguish 
them  from  the  actual  components  of  the  combined  final  rise  in  (31)  and  (32). 
Following  this  notation,  the  approximate  sum  of  cubes  of  individual  final 
rises  is  Ahf,  defined  as 


What  makes  (030)  exact  while  (033)  is  only  an  aproximation?  The  answer 
is  that  in  the  actual  combined  plume  there  is  a single  time  to  final  rise  t^ 
that  applies  to  both  the  momentum  and  buoyancy  components.  In  contrast  the 
individual  plumes  in  the  approxima tion  of  (033)  each  have  their  own  final 
rise  times  and  t^^  . The  approximation  lies  in  adding  events  that  occur 
at  two  different  times  using  a sum  of  cubes  law  that  is  really  only  valid  for 
superposition  of  effects  at  the  same  time  of  travel. 

To  find  the  ratio  the  exact  combined  rise  in  (030)  to  the  approximate 


sum  of  cubes  law  to  the  final  rises  of  the  momentum  jet  and  buoyant  plume 


Ahf  (^^Mf  ^ 


(033) 


rise  in  (033)  we  first  factor  out  the  momentum  rise  components  to  write 
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This  ratio  will  depend  on  the  way  in  which  the  atmospheric  turbulence 
dissipation  e varies  with  height.  Using  the  definitions  of  the  rise 
components  in  (D16),  (017)  and  (031)  and  (032)  leads,  after  some  algebraic 
manipulation,  to  the  final  forms 

AHf  [ 1 + I Tj  ] 

= for  e E constant  (035) 


where  ^ and  r^  are  defined  by  (023),  (024),  (025)  and  (026)  for  constant  e. 
This  constant  value  of  e i s a reasonable  approximation  for  high  stacks  or  for 
convective  turbulence. 

For  the  other  limiting  case  of  a ground  level  release  into  mechanical 
turbulence,  (034)  and  the  definitions  of  the  rise  components  lead  to 

for  E « (Ah)'^  (036) 


Ahf  0+1 
Ahf  [1+5 


1/2 


3i  1/3 


where  t and  the  time  ratio  r^  are  defined  by  (023),  (027),  (028)  and  (029) 
for  release  into  turbulence  dissipation  e that  varies  as  (Ah)"^. 
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